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Topological insulators represent a new quantum state of matter which is characterized by peculiar edge or surface 
states that show up due to a topological character of the bulk wave functions. This review presents a pedagogical account 
on topological insulator materials with an emphasis on basic theory and materials properties. After presenting a historical 
perspective and basic theories of topological insulators, it discusses all the topological insulator materials discovered to 
date, with some illustrative descriptions of the developments in materials discoveries in which the author was involved. 
A summary is given for possible ways to confirm the topological nature in a candidate material. Various synthesis 
techniques as well as the defect chemistry that became important for realizing bulk-insulating samples are discussed. 
Characteristic properties of topological insulators are discussed with an emphasis on transport properties. In particular, 
the Dirac fermion physics and the resulting peculiar quantum oscillation patterns are discussed in detail. It is emphasized 
that proper analyses of quantum oscillations make it possible to unambiguously identify surface Dirac fermions through 
transport measurements. The prospects of topological insulator materials for elucidating novel quantum phenomena that 
await discovery conclude the review. 
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1. Introduction 

The progress in condensed matter physics is often driven by 
discoveries of novel materials. In this regard, materials pre- 
senting unique quantum-mechanical properties are of particu- 
lar importance. Topological insulators (TIs) are a class of such 
materials and they are currently creating a surge of research 
activities. ^~^^ Because TIs concern a qualitatively new aspect 
of quantum mechanics, i.e. the topology of the Hilbert space, 
they opened a new window for understanding the elaborate 
workings of nature. 

TIs are called "topological" because the wave functions 
describing their electronic states span a Hilbert space that 
has a nontrivial topology. Depending on the way the Hilbert- 
space topology becomes nontrivial, there can be various dif- 
ferent kinds of TIs."^^ An important consequence of a nontriv- 
ial topology associated with the wave functions of an insulator 
is that a gapless interface state necessary shows up when the 
insulator is physically terminated and faces an ordinary insu- 
lator (including the vacuum). This is because the nontrivial 
topology is a discrete characteristic of gapped energy states, 
and as long as the energy gap remains open, the topology can- 
not change; hence, in order for the topology to change across 
the interface into a trivial one, the gap must close at the inter- 
face. Therefore, three-dimensional (3D) TIs are always asso- 
ciated with gapless surface states, and so are two-dimensional 
(2D) TIs with gapless edge states. This principle for the neces- 
sary occurrence of gapless interface states is called bulk-edge 
correspondence in topological phases. 

A large part of the unique quantum-mechanical proper- 
ties of TIs come from the peculiar characteristics of the 



edge/surface states. Currently, the TI research is focused 
mostly on time-reversal (TR) invariant systems, where the 
nontrivial topology is protected by time-reversal symmetry 
(TRS).^~^^ In those systems, the edge/surface states present 
Dirac dispersions (Fig. 1), and hence the physics of relativistic 
Dirac fermions becomes relevant. Furthermore, spin degener- 
acy is lifted in the Dirac fermions residing in the edge/surface 
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Fig. 1. (Color online) Edge and surface states of topological insulators 
with Dirac dispersions, (a) Schematic real-space picture of the ID helical 
edge state of a 2D TI. (b) Energy dispersion of the spin non-degenerate edge 
state of a 2D TI forming a ID Dirac cone, (c) Schematic real-space picture 
of the 2D helical surface state of a 3D TI. (d) Energy dispersion of the spin 
non-degenerate surface state of a 3D TI forming a 2D Dirac cone; due to the 
helical spin polarization, back scattering from k to -k is prohibited. 
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States of TR-invariant TIs and their spin is locked to the mo- 
mentum (Fig. 1). Such a spin state is said to have "helical 
spin polarization" and it provides an opportunity to realize 
Majorana fermions^^ in the presence of proximity-induced su- 
perconductivity,^^ not to mention its obvious implications for 
spintronics applications. An important objective of the exper- 
imental studies of TIs has been to elucidate the existence and 
the nature of such helically spin-polarized Dirac fermions in 
the topological surface state. 

This review is aimed at providing a pedagogical introduc- 
tion to the field of topological insulators, putting emphasis on 
the basic theory and materials properties. It also elaborates on 
the basics of the characterizations of TI materials using trans- 
port measurements. 

2. Historical Perspective 

It is useful to understand the historical background in or- 
der to appreciate the importance of TIs in condensed matter 
physics. In this section, I will briefly discuss several relevant 
issues that preceded the discovery of TIs and also describe 
how the TI research developed in its early days. 

2. 1 Integer quantum Hall effect 

In 1980, von Klitzing et ah discovered the quantum Hall 
efl'ect in a high-mobility 2D semiconductor under high mag- 
netic fields. ^^ The occurrence of this eff'ect is usually limited 
to very low temperatures, where localization of electrons and 
Landau quantization of their energy spectrum lead to van- 
ishing longitudinal conductivity cTxx together with quantiza- 
tion of the Hall conductivity CFxy to integer multiples of e^ Ih 
when the chemical potential is located in between Landau 
levels. Such a quantization of transport coefficients obviously 
pointed to a macroscopic quantum phenomenon, and in 1982, 
it was recognized by Thouless, Kohmoto, Nightingale, and 
den Nijs (TKNN)^^ that this phenomenon not only is quantum 
mechanical but also is topological, namely, TKNN showed^^ 
that the Hilbert space spanned by the wave functions of a 
quantum Hall system obtains a nontrivial topology, which can 
be quantitatively specified by an integer topological invariant 
called TKNN invariant v, and that cr^y becomes equal to v 
times e^ Ih. The TKNN invariant is also called the first Chern 
number or the winding number, and it is equal to the Berry 
phase of the Bloch wave function calculated around the Bril- 
louin zone (BZ) boundary divided by In (actual calculations 
are shown in the next section). 

In hindsight, the quantum Hall system can be considered to 
be the first topological insulator that became known to physi- 
cists, because when the quantization is taking place, the en- 
ergy spectrum is gapped due to the Landau quantization and 
the chemical potential is located within the gap, which is a sit- 
uation akin to an insulator. In this case, the nontrivial topology 
specified by the TKNN invariant is characteristic of a 2D sys- 
tem with broken TRS. Also, as was shown by Halperin,^^ the 
integer quantum Hall eff'ect is always accompanied by chiral 
edge states, and those gapless states residing at the interface to 



the vacuum can be understood to be a result of the bulk-edge 
correspondence due to the topological 2D "bulk" state. 

2.2 Quantum spin Hall effect and Z^ topology 

On a diff'erent front in condensed matter physics, genera- 
tion and manipulation of spin currents have been attracting 
a lot of interest, since they will have a profound impact on 
future spintronics. ^^^ In this regard, the spin Hall eff'ect, the 
appearance of transverse spin current in response to longi- 
tudinal electric field, has been discussed theoretically since 
1970s, ^^~^^^ but its experimental confirmation by Kato et al. ^^^ 
in 2004 gave a big boost to the research of this phenomenon. It 
was soon recognized that the spin Hall eff'ect in nonmagnetic 
systems is fundamentally related to the anomalous Hall eff'ect 
in ferromagnets,^^^ and similarly to the latter eff'ect, there are 
both intrinsic and extrinsic origins of the spin Hall eff'ect. The 
intrinsic mechanism of the spin Hall eff'ect stems from the 
Berry curvature of the valence-band Bloch wave functions in- 
tegrated over the Brillouin zone.^^'^^^ Since such an integral 
can become finite even in an insulator, Murakami, Nagaosa, 
and Zhang went on to propose the idea of spin Hall insula- 
tor }'^^ which is a gapped insulator with zero charge conduc- 
tivity but has a finite spin Hall conductivity due to a finite 
Berry phase of the occupied states. 

Later it was shown^^^ that the proposed spin Hall insulators 
cannot really generate spin currents in the absence of any elec- 
trons at the Fermi level, but this idea triggered subsequent pro- 
posals of its quantized version, the quantum spin Hall (QSH) 
insulator, by Kane and Mele,^^'^^^ followed by an independent 
proposal by Bernevig and Zhang. ^^^ The QSH insulators are 
essentially two copies of the quantum Hall system, in which 
the chiral edge state is spin polarized and the two states form 
a time-reversed pair to recover the overall TRS. When current 
flows using the edge states of a QSH insulator, a quantized 
version of the spin Hall efl'ect, QSH efl'ect, is predicted to be 
observed. Since the predicted phenomenon is based on the 
quantum Hall efl'ect, it only exists in 2D. While it is not a 
priori clear how one can achieve such a state with quantized 
edge states in zero magnetic field, the ingenious proposal by 
Kane and Mele provided a concrete model to realize the QSH 
insulator. ^^^ Their model is essentially a graphene model with 
spin-orbit coupling (SOC). 

In graphene, the band structure near the Fermi level con- 
sists of two linearly dispersing cones located at K and K' 
points in the BZ;^^^ since the low-energy physics on these 
cones is described by employing the Dirac equation with 
the rest mass set to zero,^"^^ this dispersion is called Dirac 
cone and the electrons are said to behave as massless Dirac 
fermions. Kane and Mele showed^^^ that a finite SOC leads to 
an opening of the gap at the crossing point of the cone (called 
Dirac point) and, furthermore, that a time-reversed pair of 
spin-polarized one-dimensional (ID) states indeed show up at 
the edge in some parameter range; in this model, the desired 
spin polarization of the edge state is achieved due to the SOC 



J. Phys. Soc. Jpn. 



INVITED REVIEW PAPERS 



which has an inherent tendency to aUgn spins in relation to 
the momentum direction. This pecuHar spin-non-degenerate 
state [Fig. l(a,b)] is often said to have helical spin polar- 
ization or spin-momentum locking. Intriguingly, those elec- 
trons in the gapless edge state behave as ID massless Dirac 
fermions within the gap opened in the 2D Dirac cone. In this 
case, the 2D "bulk" electrons can be viewed as massive Dirac 
fermions because of the finite energy gap at the Dirac point. 

Most importantly, Kane and Mele recognized that the elec- 
tronic states of their QSH insulator is characterized by a novel 
topology specified by a Z2 index,^^^ which expresses whether 
the number of times the ID edge state crosses the Fermi level 
between and n/a is even or odd (a is the lattice constant). 
Remember, in mathematics the group of integer numbers is 
called Z and its quotient group classifying even and odd num- 
bers is called Z2 ; hence, a Z2 index generally gives a topolog- 
ical classification based on parity. (A detailed description of 
the Z2 index for TR-invariant TIs is given in the next section.) 
The theoretical discovery of the Z2 topology in insulators was 
a big step in our understanding of topological phases of mat- 
ter, because it indicated that nontrivial topology can be em- 
bedded in the band structure of an ordinary insulator and that 
breaking of TRS is not mandatory for realizing a topological 
phase. 

Unfortunately, the SOC in graphene is very weak, and 
hence it is difficult to experimentally observe the QSH eff'ect 
predicted in the Kane-Mele model. However, another theoret- 
ical breakthrough was soon made by Bemevig, Hughes, and 
Zhang (BHZ),^^^ who constructed a 2D model to produce a 
Z2 topological phase based on the band structure of HgTe; 
based on their model, BHZ predicted that a CdTe/HgTe/CdTe 
quantum well should give rise to the QSH eff'ect. This pre- 
diction was verified in 2007 by Konig et al. ,^^^ who observed 
CTjcx to be quantized to 2e^/h in zero magnetic field when the 
chemical potential is tuned into the bulk band gap, giving evi- 
dence for the gapless edge states in the insulting regime. This 
was the first experimental confirmation of the TR-invariant TI 
characterized by the Z2 topology. 

Without waiting for experimental verification of the Z2 
topology in 2D, theorists noticed that this topological classi- 
fication of insulators can be extended to 3D systems, where 
there are four Z2 invariants to fully characterize the topol- 
ogy.^^~^^^ In fact, the term "topological insulator" was coined 
by Moore and Balents in their paper to propose the exis- 
tence of TIs in 3D systems. ^^^ For 3D TIs, Fu and Kane 
made a concrete prediction in 2006 that the Bii_;cSb;c alloy 
in the insulating composition should be a TI, and they fur- 
ther proposed that the nontrivial topology can be verified by 
looking at the surface states using the angle-resolved photoe- 
mission spectroscopy (ARPES) and counting the number of 
times the surface states cross the Fermi energy between two 
TR-invariant momenta. ^^^ The proposed experiment was con- 
ducted by Hsieh et al who reported in 2008 that Bii_;cSb;c is 
indeed a 3D TI.^^^ The experimental identification of Bii-^^Sbjc 



as a TR-invariant TI opened a lot of new experimental oppor- 
tunities to address a topological phase of matter. For exam- 
ple, the first transport study of Bii_;cSb;c to detect topological 
2D transport channels was reported by Taskin and Ando,^^^ 
and the first scanning tunneling spectroscopy (STS) study 
that addressed the peculiar spin polarization was reported by 
Roushan et al.P'' both in 2009. Direct observation of the he- 
lical spin polarization of the surface states in Bii-^^Sbj^ us- 
ing spin-resolved ARPES was first partially done by Hsieh et 
al?^"^ and then fully accomplished by Nishide et al. (Fig. 2).^^^ 

2.3 Topological field theory 

While the notion of topological insulators became popu- 
lar only after the discovery of the Z2 topology by Kane and 
Mele,^^^ there had been theoretical eff'orts to conceive topo- 
logical states of matter beyond the scope of the quantum Hall 
system. In this respect, an important development was made 
in 2001 by Zhan and Hu, who generalized the 2D quantum 
Hall state to a four-dimensional (4D) TR-invariant state pos- 
sessing an integer topological invariant. ^^^ The eff'ective field 
theory for this 4D topological system was constructed by 
Bernevig et al.^^^ After the Z2 topology was discovered for 
TR-invariant systems in 2D and 3D,^^'^^~^^^ it was shown by 
Qi, Hughes, and Zhang^^^ that the framework of topological 
field theory is useful for describing those systems as well, and 
they further demonstrated that the Z2 TIs in 2D and 3D can 
actually be deduced from the 4D eff'ective field theory by us- 
ing the dimensional reduction. From a practical point of view, 
the topological field theory is suitable for describing the elec- 
tromagnetic response of TIs and has been used for predicting 
novel topological magnetoelectric eff'ects.^^"^^^ 
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Fig. 2. (Color online) Spin-resolved surface band structure of Bii-j^Sbjc 
(x = 0.12 - 0.13) on the (111) cleaved surface. Its surface BZ and the mo- 
mentum direction of the data are shown in the inset. Dispersions shown by 
symbols are determined from the peak positions of the spin-resolved energy 
dispersion curves and they are plotted on top of the spin-integrated ARPES 
data shown in gray scale. There are three Fermi-level crossings of the surface 
states, which signifies the Z2 -nontrivial nature of this system. Taken from 
Ref. 35; copyright American Physical Society (2009). 
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2.4 Dirac materials 

It is interesting to note that the first 3D TI material 
Bii_;cSb;c. whose topological surface state consists of 2D 
massless Dirac fermions, has long been known to possess pe- 
culiar band structure to give rise to 3D massive Dirac fermions 
in the bulk."^^^ This situation is similar to the Kane-Mele 
model where ID massless Dirac fermions emerges out of 2D 
massive Dirac fermions. Since Dirac fermions play important 
roles in TIs, it would be useful to mention the history of Dirac 
physics in condensed matter. 

The semimetal Bi has been an important testing ground of 
quantum physics ever since the Shubnikov-de Haas (SdH) and 
de Haas-van Alphen (dHvA) oscillations were discovered in 
Bi in 1930.^^^ This is essentially because the extremely low 
carrier density (-10"^ per atom) and the very long mean free 
path (reaching ~1 mm) easily put the system in the "quantum 
limit" at relatively low magnetic fields. "^"^^ In the mid-20th cen- 
tury, one of the long-standing puzzles in Bi was its large dia- 
magnetism, which defies the common wisdom for magnetism 
in metals involving Pauli paramagnetism and Landau diamag- 
netism."^^^ Intriguingly, in Bii-j^Sbj^ at low Sb concentration, 
the carrier density becomes even lower than in Bi, while the 
diamagnetic susceptibility increases, which is also opposite 
to the expectation from Landau diamagnetism."^^^ To address 
the unusual electronic properties of Bi, an eff'ective two-band 
model was constructed by Cohen and Blount in 1960.^^^ In 
1964, Wolf recognized that this two-band model can be trans- 
formed into the four-component massive Dirac Hamiltonian, 
and he presented an elegant analysis of the selection rules us- 
ing the Dirac theory.^^^ This was the beginning of the notion 
of Dirac fermions in solid states. 

An important aspect of the Dirac physics is that the "inter- 
band eff'ect" (i.e. interactions between upper and lower Dirac 
cones) becomes essential. In fact, the Dirac formalism allows 
one to naturally include the interband eff'ect into calculations, 
and in 1970 Fukuyama and Kubo showed that the large dia- 
magnetism in Bi and Bii_;cSb;c is a consequence of their Dirac 
nature."^^^ 

For graphene, the mapping of the k • p Hamiltonian"^^^ to the 
massless Dirac Hamiltonian was used by Semenoff'in 1984,^^^ 
and now with the experimental realization of graphene,^^^ this 
system has become a prototypical Dirac material. One of the 
distinguishing properties of massless Dirac fermions is the 
Berry phase of tt; an important consequence of the n Berry 
phase in the condensed matter setting is the absence of back 
scattering, which was pointed out first by Ando, Nakanishi, 
and Saito in 1998.^9) 

2. 5 Old proposals for interface Dirac fermions 

It is prudent to mention that already in 1985, the appear- 
ance of 2D massless Dirac fermions at a certain interface 
was anticipated.^^' ^^^ Specifically, Volkov and Pankratov dis- 
cussed that the relative band inversion occurring in SnTe/PbTe 
and in CdTe/HgTe would lead to a gapless interface state. ^^^ 



Also, Fradkin, Dagotto, and Boyanovsky discussed in 1986 
that an antiphase domain wall in PbTe (which creates the 
atomic arrangements of Pb-Te-Pb-Pb-Te-Pb or Te-Pb-Te-Te- 
Pb-Te) can be viewed as a parity anomaly in 3 -hi dimension 
and would lead to a gapless interface state. ^^^ 

What distinguishes the recent discussions of the topolog- 
ical surface states from those old works is the notion of 
topological protection, which guarantees the stability of the 
gapless state. For example, SnTe is a band-inverted semi- 
conductor, but its band inversion takes place at even num- 
ber of points in the Brillouin zone, which makes 3D SnTe 
topologically-trivial in the Z2 topological classification. ^^^ 
This means that, even if the interface state is expected to 
appear at the SnTe/PbTe interface as a result of band inver- 
sion, its gapless nature is not protected by the Z2 topology. 
(Nevertheless, SnTe was recently found to be a topological 
crystalline insulator, ^^'^"^^ and its gapless surface state is actu- 
ally topologically protected by mirror symmetry, as discussed 
later.) In the case of the antiphase domain wall in PbTe, de- 
tailed symmetry arguments highlighted the difficulty in real- 
izing the parity anomaly.^^^ 

3. Basics of the Topology in Topological Insulators 

3. 1 Berry phase 

The notion of Berry phase^^^ is important in the discussions 
of topological phases. Here, I briefly discuss its definition and 
meaning. 

Let R(0 be a set of time-dependent parameters and consider 
it as a vector in the parameter space. We write the Hamiltonian 
of a system specified by the parameter R(0 as //[R(0] and its 
n\h eigenstate as \n, R(0). The Schrodinger equation for this 
system is 

if [R(0] \n, R(0) = ^«[R(0] \n, R(0>. (1) 

Suppose that R changes adiabatically from the ^ = value Rq. 
The time evolution of the state follows 



H[^{t)]\nj) = in-\nj) 
ot 



(2) 



and the state at time t is given by 

\n, t) = exp U J dfLniW)]] \n, R(0>, (3) 

where 

LJR(0] = inR(t) . {n, R(t)\WR\n, R(0) - En[R(t)l (4) 

This can be easily confirmed by putting Eq. (3) into the r.h.s. 
of Eq. (2). Using Ln[R(t)] given in Eq. (4), one may write the 
time-dependent state as 

\n, t) = exp I - r dfR(f) • {n, R(f)\VR\n, R(f)}\ \n, R(0) 

xexpl^ I dfEn[R(f 



)] . (5) 
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In this expression of \n, t), the first exponential term repre- 
sents the nontrivial eff'ect of the quantum-mechanical phase 
accumulated during the time evolution, and the last exponen- 
tial term is a trivial one called dynamical term. 

When R moves on a closed loop C from ^ = and returns 
to the original position at t = T, i.e. R(T) = Rq, the Berry 
phase YniC] for this loop C is defined as 

yAC] ^ I J^R(0-/<^,R(OIV/^KR(0> (6) 



-i 



JR-/<n,R|V/?KR) 



i- 



= -(hdR- A„(R) 



■X 



dS ■ B„(R) 



(7) 



(8) 



(9) 



The last equality comes from the Stokes' theorem. Here, we 
define the Berry connection 



A„(R) = -/<^,R|V^KR), 
and its rotation is the Berry curvature 

B„(R) = V^xA„(R). 



(10) 



(11) 



From Eqs. (5) and (6), one can see that the Berry phase means 
the accumulated phase factor of a quantum-mechanical sys- 
tem after it completes a closed path in the parameter space. 
The Berry connection corresponds to the gauge field defined 
on that parameter space, similar to the vector potential for 
electromagnetic fields in real space. 

3.2 TKNN invariant 

The topological invariant defined for the integer quantum 
Hall system, the TKNN invariant,^^ is closely related to the 
Berry phase. To see this, we derive the TKNN invariant by 
calculating the Hall conductivity of a 2D electron system of 
size L X L in perpendicular magnetic fields, where the electric 
field E and the magnetic field B are applied along the y and 
z axes, respectively. By treating the eff'ect of the electric field 
£^ as a perturbation potential V = -eEy, one may use the 
perturbation theory to approximate the perturbed eigenstate 
\n)E as 



\n)E = 1^) + ^ 



{m\ j-eEy) \n) 



m) + 



(12) 



Using this perturbed eigenstate, one may obtain the expecta- 
tion value of the current density along the x axis, jx, in the 
presence of the E field along the y axis as 



{Jx)E = Yjf^^-^^''\^ijj)\''^^ 



{Jx}e^O + 






{n\ evx \m}{m\ (-eEy) \n} 



{n\ (-eEy) \m}{m\ eVx \n) 



(13) 



where Vx is the electron velocity along the x direction and 
f{En) is the Fermi distribution function. The Heisenberg 
equation of motion j^y = Vy = ^ [y, H] leads to 



{m\ Vy \n) = —(En- Ejn){m\y\n), 



''- in 

which allows one to evaluate 

{Jx)e 



(14) 



E 

ifie^ 
1? 



Yjf(^n) 



{n\ Vx \m){m\ Vy \n} - {n\ Vy \m}{m\ Vx \n} 
X J- . (15) 

When we consider a system in a periodic potential and its 
Bloch states |w„k) as the eigenstates, the identity 

1 d 

in okj^ 

allows one to rewrite Eq. (15) into the form 



cr 



xy 



L2 



k n^m 



id d d d \ 

Since the Berry connection defined in Eq. (10) is written for 
Bloch states as 

a„(k) = -i{Un^\ "Vk \Unk) = -i{Unk\ -^ \Unk), (18) 

dk 
the Hall conductivity reduces to 



(Trv = V — 



with 






(19) 



(20) 



This y can be expressed as v = 2n ^n with Vn the contribution 
from the ^th band, and one can easily see that y„ is related to 
the Berry phase, namely 

r d^k i danj _ dam A 
Jbz 271 \ dkx dky ) 

= —d) dk- a„(k) 

^^ JdBZ 



= -lyJ^BZ]. 

2n 



(21) 



Because of the single- valued nature of the wave function, its 
change in the phase factor after encircling the Brillouin zone 
boundary (^BZ) can only be an integer multiple of In, which 
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Fig. 3. (Color online) Kramers pairs of bands. Note that each pair of bands 
are degenerate at the TR-invariant momentum (TRIM) where +k becomes 
equivalent to -k due to the periodicity of the Brillouin zone. In this figure, 
there are two TRIMs, k = and k = n (which is equivalent to /: = -n). The 
lifting of degeneracy at k values other than and n comes from spin-orbit 
coupling. 



means 



ynWBZ] = 2nm (m e Z). 



(22) 



Therefore, y„ can only take an integer value, and hence CTxy is 
quantized to integer multiples of e^/h. The integer y is called 
TKNN invariant, and it plays the role of the topological in- 
variant of the quantum Hall system, which is a TRS -breaking 
TI. 

3.3 Time-reversal operator 

We now consider insulators preserving TRS. The TR oper- 
ator for a spin 1/2 particle takes the simple form = -iSyK, 
where K is the complex conjugate operator. Hereafter, Sj^ 
(ji = x,y, z) denotes the spin operator given by Pauli matri- 
ces. An important property of the TR operator is that 

^2 



0" 



-1. 



(23) 



Taking the eigenstates of s^ as the basis set { |(t) }, one may 
confirm 

<^|0|0) = - YjWcT}{cT\iSy\cT'}{cT'\cP*} 



= -<0l©l^>. 

Similarly, it is easy to confirm 

and for arbitrary linear operator A 

{@iJ/\@A@-^\@4>} = {4>\A^ifr}. 



(24) 



(25) 



(26) 



3.4 TRS and Block Hamiltonian 

Let "H be the total Hamiltonian of a periodic system: 

'H\^n\d = En^,\l|fn^^• (27) 

B loch's theorem leads to the separation of </^„k into 

Wny.) = ^'"^"l^nk), (28) 
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Fig. 4. (Color online) Time-reversal-invariant momenta (TRIMs). (a) 
There are 4 TRIMs for a 2D BZ; this figure shows the case of a square BZ. (b) 
The TRIMs for a TR-symmetric ID system considered in Sec. 3.5. One can 
see that identifying t ^ ky maps the periodic (k, t) space to the ordinary 2D 



BZ shown in (a), (c) There are l 
the case of a cubic BZ. 



TRIMs for a 3D BZ; this figure exemplifies 



where |w„k) is the cell-periodic eigenstate of the Bloch Hamil- 
tonian 



/kr (Tj J\iY 



H{\i) = e-'^' "H e' 
and |w„k) satisfies the reduced Schrodinger equation 

if(k)|w„k) = En^\Un^). 



(29) 



(30) 



It is important to note that when "H preserves TRS, i.e. 
['K, 0] = 0, then if (k) satisfies 

if(-k) = 07f(k)0-^ (31) 

This identity means that the energy bands of a TR sym- 
metric system come in pairs (i.e. -hk state and -k state are at 
the same energy), and they are called Kramers pairs (Fig. 3). 
Naturally, Kramers pairs are degenerate at the TR-invariant 
momentum (TRIM) where -hk becomes equivalent to -k due 
to the periodicity of the BZ (Fig. 4). 

Now we discuss a matrix representation of the TR operator 
in the Bloch wave function basis. A convenient matrix is 



WQ,^(k) = <w^,-k|0|W;S,k), 



(32) 



where a and p are band indices. This matrix relates the two 
Bloch states \ua-\i) and |w^,k) via 



\Ua-^) = ^ Kfi(^) ® I^Ak>- 



(33) 
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One can easily confirm that Wafi(k) is a unitary matrix, i.e. 

^ wl^(k)wa^(k) = <w^,k|w^,k) = ^^y, (34) 

a 

and it has the following property: 

w^^(-k) = -w^^(k). (35) 

This equation means that at a TRIM A, the w matrix becomes 
antisymmetric, i.e. 

^/3aW = -Wa/sW. (36) 

For example, if there are only two occupied bands, Wa^ be- 
comes a 2 X 2 matrix and at A it can be explicitly written 
as 



w(A): 



wi2(A) \ ,.J 1 



(37) 

Another convenient matrix to consider for a TR symmetric 
system is the U(2) Berry connection matrix 

aa^(k) = -/<Wa,k|Vk|w^,k). (38) 

Notice that a is actually a set of three matrices. By using Eqs. 
(25) and (33), it can be shown that Sia^ik) and a^^i-k) are 
related by 

a(-k) = w(k)a*(k)w^(k) + /w(k)Vkw"^(k). (39) 

By taking the trace and replacing -k -^ k, one obtains 

tr[a(k)] = tr[a(-k)] + /tr[w^(k)Vkw(k)]. (40) 

This relation becomes important in the calculation of the Z2 
topological invariant. 

3.5 Z2 time-reversal polarization 

In order to pedagogically derive the topological invariant 
to characterize 2D electron systems preserving TRS, we fol- 
low Fu and Kane^^^ and consider a ID system with length 
L and lattice constant a =1. We only consider two bands that 
form a Kramers pair, and denote their Bloch wave functions as 
\ui(k)} and \u2(k)}. Assume that the band parameters change 
with time and return to the original values at ^ = T. Further- 
more, we consider the situation when the ID Hamiltonian sat- 
isfies the following conditions: 



H[t -\-T]= H[T] 
H[-t] = @H[t]@-\ 



(41) 
(42) 



From the modem theory of ferroelectricity,^^^ it is known 
that charge polarization Pp can be calculated by integrating 
the Berry connection of the occupied states over the BZ. In the 
present case of the ID two-band system, Pp may be written as 

dk 



-f 



In 



Aik\ 



(43) 



where 

A{k) = -i{ui(k)\Vk\ui(k)}-i{u2(k)\Vk\u2(k)} (44) 

= an(k) -\- a22(k) 

= trM. (45) 

The contribution from each band may be called partial polar- 
ization defined by 



r dk 



(k) (/=1,2), 



(46) 



which gives the relation Pp = Pi -\- P2. The time-reversal po- 
larization is defined by 



Pe = Pi-P2 = 2P,-Pp. 



(47) 



Intuitively, Pq gives the diff'erence in charge polarization be- 
tween spin-up and spin-down bands, since \ui(k)} and \u2(k)} 
form a Kramers pair. 

From Eqs. (41) and (42), one can see that the system is TR 
symmetric at ^ = and t = T/2. At these times, the Kramers 
degeneracy must be observed at any k, which dictates that 
the time-reversed version of \u2(k)), i.e. 0|w2(^)), is equal to 
\ui(-k)} except for a phase factor. Hence, at ^ = and t = T/2, 

e\u2(k)} = e-'^^%i(-k)} 

&\ui(k)) = -e-'^^-%2(-k)). 



(48) 
(49) 



Using these relations, the w matrix can be shown to become 



w(k) = 







-ix(k) 



-e-^^i-k) 



(50) 



Now we calculate Pi at the TR symmetric times. First, us- 
ing Eq. (49) one may obtain 



which leads to 
Pi = 



an(-k) = a22(k) - -^x(k), 
ok 



— \ I dkau(k)-\- I dka 

'^ \ Jo J-n 

1 r / 

—I dk\a 



(51) 



i(k) 



nik) + a22ik) - ^X(k) 



A(k)-—{x{7i)-x{0)l (52) 

Jo 27r 271 

Since wuik) = ^"'^^^^ from Eq. (50),;^(^) can be written as 

X(k) = ilogwnik) (53) 

and Eq. (52) reduces to 

Jo 27T 271 Wi2(0)' 

This expression leads to 

Pe = 2Pi-Pp= r ^[A(k)-A(-k)]--log^^^. (55) 
Jo 2;: 71 wi2(0) 
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Noticing A(k) = tr[a(^)] and using Eq. (40), one obtains 
— itr[w\k)—-w(k)] - - log — 

ZTT ok 71 Wi2(0) 

— — log(det[wW]) - - log -— 

ZTT ok 71 Wi2(0) 



- /" 

Jo 

- '/ 

Jo 
i 

71 

-..,2 



1 ^ det[w(7r)] 
- log 

2 ^ det[w(0)] 



/ wi2(7r) 
- log . 

71 Wi2(0) 



(56) 



Since det[w] = w\^ in the present case, Eq. (56) reduces to 



P,^ 



1 



l0£ 



m 



Vwi2(0)2 Wi2(7r) 



>^12(0) V>^12W^J 



(57) 



The argument of this log is +1 or -1. Since log(-l) = i7i if 
the angle in the complex plane is restricted from to l7i, one 
can see that Pq is or 1 (mod 2). Physically, the two values of 
Pq corresponds to two different polarization states which the 
present ID system can take at ^ = and t - T 12. 

We now consider the change of Pq from ^ = to T 12. The 
wave function \unik, t)) can be viewed as a map from the 2D 
phase space (k, t), which forms a torus due to the periodic 
boundary conditions, to the Hilbert space. This Hilbert space 
can be classified into two groups depending on the diff'erence 
in Pq between ^ = and t = T/2, 



A = Pe(T/2)-Pe(0). 



(58) 



This A is specified only in mod 2, so it gives a Z2 topolog- 
ical invariant to characterize the Hilbert space. Intuitively, 
when Pq changes between t = to T/2, the Hilbert space is 
"twisted" and A = 1, while the Hilbert space is trivial (A = 0) 
when there is no change in Pq. Using Eq. (57), A can be given 
in terms of 



Wl2(A/) 



(-i)^=n^^^^^' 



(59) 



where Ai = (k, t) = (0, 0), A2 = (71, 0), A3 = (0, T/2), and A4 = 
(71, T/2) [see Fig. 4(b)]. The physical consequence of a cycle 
with A = 1 is spin pumping from one end of the ID system to 
the other.^^^ 

3.6 General formula of the Z2 invariant 

Extension of the above argument to a multiband system is 
not very difficult. The Hamiltonian still satisfies Eqs. (41) and 
(42), and in the following we take T = 27: for simplicity. 
Consider that 2N bands are occupied and forming N Kramers 
pairs. For each Kramers pair n, at the TR symmetric times 
^ = and 71 the wave functions are related by 

0|w^(yl)) = e-'^'^^%1(-k)) 



^M(k)} 



-iXn(-k)\j.n 



4(-k)), 



(60) 
(61) 








g-'Xi(k) 








iXll 


-k) 






















g-iX2(k) 










-g->X2(-k) 






and the w matrix is given by 



w(k) = 



(62) 
Hence, at ^ = and 71, w(0) and w(7i) become antisymmetric 
and we have 

W12(A0W34(A0 . . . W2A^-l,2A^(A0 = ^-^"^^i^-^^') 

= Pf[w(A,)]. (63) 

Note that in the above equation, w is viewed as a function of 
k and t, and we used the formula for the Pfaffian of a 2N x 2N 
skew- symmetric matrix with 2x2 blocks on the diagonal. In 
general, Pfaffian is defined for an antisymmetric matrix and is 
related to the determinant by 



Ff[Af = det[A]. 



(64) 



It is straightforward to extend the previous calculations for 
the TR symmetric times ^ = and 71 to obtain 



^. = rs^w-sB^* 



)-xnm 



r dk i / Pf[w(7r)] 

which in turn gives the TR polarization 



(65) 



^^ / Vdet[w(Q)] Pf[w(7r)] 
' i7i ""^X Pf[w(0)] ' Vdet[w(7r)]'' ^ 



Therefore, the Z2 topological invariant v is given by 

' Pf[w(A,)] 



(-1)^ = n 



\J^ Vdet[w(A,)]' 



(67) 



which classifies the Hilbert space into "twisted" (v = 1) and 
trivial (v = 0) ones. 

By reinterpreting the periodic 2D phase space {k, t), which 
forms a torus, as the 2D Brillouin zone (kx,ky), the above 
theory provides a Z2 topological classification of 2D TR- 
invariant insulators with 2N occupied bands. 

3. 7 Extension to 3D systems 

Based on a general homotopy argument, Moore and Ba- 
lents showed^^^ that there are four Z2 invariants for 3D sys- 
tems. While the construction of the homotopy was involved, 
the physical origin of the four invariants can be easily under- 
stood. For simplicity, consider a cubic system and take the 
lattice constant a = 1. In the 3D BZ of this system, there 
are eight TRIMs denoted as Ao,o,o, K,o,o, ^o,7t,o, Ao,o,;r, K,o,7:, 
Ao,;r,;r, K,7:,o, and A;r,;r,;r [sec Fig. 4(c)]. At thcsc points in 
the BZ, the Bloch Hamiltonian becomes TR symmetric, i.e. 
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&H(Ai)&~^ = H(Ai), and each Kramers pair of bands be- 
come degenerate. 

Notice that the six planes in the 3D BZ, x = 0, x = ±7i, 
y = 0, y = ±71, z = 0, and z = ±7i possess the symmetries 
of the 2D BZ, and therefore they each have a Z2 invariant. 
The six invariants may be denoted as xq, xi, yo, yi, zo, and 
Zi, but those six invariants are not all independent. ^^^ This 
is because the products xqXi, yoyi, and zoZi are redundant, 
which stems from the fact that those three products involve 
Pf[w(A/)]/ Vdet[w(A/)] from all eight TRIMs and hence are 
the same. This means that there are two constraining relations 
•^0-^1 = yoyi = ZoZi, which dictates that only four invariants 
can be independently determined in a 3D system. 

The concrete construction of the four Z2 invariants were 
given by Fu and Kane.^^^ For each TRIM A/, we define 

Pf[w(A,)] 



S(Ai) = 



Vdet[w(A,)] ' 



(68) 



Using this (5(A,), the four Z2 invariants vo, vi, V2, V3 are de- 
fined as 



(-ir = n 5(A„„„„„3) 



(69) 



nj=0,7i 



(-IT 



n S(A„,,„„„,) (/= 1,2,3). (70) 



nj^i=0,7T;ni-7T 



The invariant vq is a product of all eight ^(A/)'s, so it is unique 
to a 3D system. On the other hand, other v/'s are a product of 
four ^(A/)'s and is similar to the Z2 invariant in the 2D case. 
For example. 



VS = ^(Ao,0,;r) ^(K,0,n) ^(Ao,;r,;r) ^(K,7:,7:) 



(71) 



corresponds to the Z2 invariant on the z = tt plane, which can 
be seen by considering the TR polarization defined on this 
plane, Pgiy = 0)^.;, = ^(Ao,o,;r) ^(A;r,o,;r) and Pe(y = n\^n = 
^(Ao,;r,;r) ^(A;r,;r,;r)- When the TR polarization changes between 
y = and y = n [i.e. Pgiy = 0)^^;^ and Pgiy = 7i)z^j^ are diff'er- 
ent], then the Z2 topology is nontrivial and vs becomes 1. 

The physical consequence of a nontrivial Z2 invariant is the 
appearance of topologically -protected surface states. This is 
graphically shown in Fig. 5, in which topologically trivial and 
nontrivial surface states are compared. In the nontrivial case, 
the Kramers pairs in the surface state "switch partners", and 
as a result, the surface state is guaranteed to cross any Fermi 
energy inside the bulk gap. This switch-partner characteristic 
is a reflection of the change in the TR polarization discussed 
above. 

It is customary to write the combination of the four invari- 
ants in the form (vq; vi V2 V3), because (vi V2 V3) can be inter- 
preted as Miller indices to specify the direction of vector A/ 
in the reciprocal space. A 3D TI is called "strong" when vq = 
1, while it is called "weak" when vq = and v/ = 1 for some / 
= 1,2,3. 



(a) 



EA 



r?s2^ 





Fig. 5. (Color online) Schematic pictures of the surface states between two 
surface TRIMs, A^ and A^, for (a) topologically trivial and (b) topologically 
nontrivial cases. In the latter, the "switch-partner" characteristic is a reflection 
of the change in the TR polarization (see text). The shaded regions represent 
the bulk continuum states. To enhance the visibility of the Kramers degener- 
acy, surface state dispersions are shown from -A^ to A^ through A^. 



3.8 Topological insulators with inversion symmetry 
Define the space inversion operator IT as 



n|x,cr) = |-X,0-), 



(72) 



where x is the spatial coordinate and a signifies the spin. One 
can easily see that 

tf = 1, (73) 

which means that the eigenvalue of the IT operator is ± 1 . It is 
easy to see that 11 operates on the wave function in momentum 
representation as n|k, cr) = |-k, cr), and hence Hamiltonian 

H = 2k,o-,o-' |k, (T)ifo-,o-'(k)(k, cr'l transforms as 

UHU-' = Yj \-K(t}H^,A\^){-K(t'\ 

k,cr,cr' 

= Yj |k,<r)7/^,<.'(-k)<k,(r'|. (74) 



k,cr,cr' 

When H satisfies H(-k) = //(k), it also satisfies 



(75) 



and hence the system preserves inversion symmetry (IS). Our 
aim is to express the Z2 invariant using the eigenvalues of IT, 
namely, the parity. 

We have previously defined the U(2) Berry connection ma- 
trix Siafiik) in Eq. (38) and calculated its trace, tr[a(k)], in Eq. 
(40). Now we denote this trace as a^(k) and define 



F(k) = V,xa^(k), 



(76) 



which is the Berry curvature of a^. The two symmetries dic- 
tate that this F(k) has the following properties :^^^ 



TRS 
IS 



F(-k) = -F(k) 
F(-k) = ±F(k). 



(77) 
(78) 



Therefore, the Berry curvature F(k) = for any k, and hence 



J. Phys. Soc. Jpn. 



INVITED REVIEW PAPERS 



one can always chose a gauge which makes a^(k) = 0. 

To actually obtain a gauge to set a^(k) = for any k, we 
consider the matrix 



VQr^(k) = (w^k|n0|t/^k). 



(79) 



One can see that this v-matrix is antisymmetric (Va^ = -v^a) 
and unitary (2^ ^afi^lg = ^a^)- Furthermore, v is related to 
a^(k) via 



^tr[v^V^v] = a^(k), 



(80) 



which can be verified by calculations similar to those em- 
ployed for deriving Eq. (40). Since tr[v^Vy^v] = Vy^trQogv] = 
V;^log(det[v]),^^^ one can see 

a^(k) = ^ log(det[v(k)]) = Nk log Pf [v(k)]. (81) 

To set a^(k) = 0, all we have to do is to adjust the phase of 
\uak} to make Pf[v(k)] = 1. 

Now we calculate the w-matrix of an inversion symmetric 
system. Let ^^(A/) = ±1 be the eigenvalue of IT for band a at 
TRIM A/, namely, Il\ua(Ai)) = ^q,(A/)|wq,(A/)). Remembering 
n^ = 1 , one may calculate 

w^^(Ai) = {u^(-Ai)\ n n K(Ai)} 

= U^i)vap{^i\ (82) 

Since w^p and v^p are antisymmetric, ^^ - ^j3 when w^js = 
-Wjsa ^ 0. Such a non-zero Wap is obtained only when the 
bands a andyS form a Kramers pair. Hence, if the bands a and 
P are the n\h Kramers pair in the total of 2N bands, we may 
write ^a = ^p = ^2n- Remembering that Wa^ has the peculiar 
skew- symmetric from of Eq. (62), one can see 

Pf [w(A,)] = Pf [v(A,)] Y] ^2n(At). (83) 

Since we have chosen the gauge to give Pf[v(k)] = 1, we fi- 
nally obtain 



S(Ai) 



Pf[w(AO] 
Vdet[w(A,)] 



n^2«(Ao. 



(84) 



n^l 



This means that the Z2 invariant can be calculated simply by 
the parity eigenvalues ^2n at TRIMs A/. For example, for a 2D 
system with 2N bands preserving IS and TRS, the Z2 invariant 
y is given by 



4 A^ 



(-1)- = P[^(A,) = P[]^^2.(A,). 



(85) 



= 1 n=l 



This simplification of the calculation of the Z2 index, accom- 
plished by Fu and Kane,^^^ greatly facilitated the theoretical 
predictions of candidate TI materials based on ab initio band 
calculations. 



3.9 BHZ model 

Now that we have learned how the topology of a TR- 
invariant insulator can be specified, let us consider a con- 
crete example that gives rise to a topologically nontrivial state. 
As mentioned in Section 2.2, the first theoretical model pro- 
posed in this context was the Kane-Mele model based on 
graphene,^^^ but it was the BHZ model^^^ that actually lead 
to the experimental discovery of the first 2D TI material. In 
this section, we discuss the BHZ model within the framework 
developed by Fu and Kane^^^ to make its Z2 topology trans- 
parent. 

The BHZ model is motivated by CdTe/HgTe/CdTe quan- 
tum well, where the s- and /^-orbital bands invert at k = as 
a function of the well width. In this quantum-well structure, 
two ^--orbital states and two /7^^3/2-orbital states appear near 
the Fermi energy. Hence, we consider the following 4 bands: 

|^,T>, \sd}, IPx'^iPyA}, \Px-iPyA-} 

Since bulk HgTe crystalizes in the cubic zinc-blende struc- 
ture, we consider 2D square lattice with two orthogonal lattice 
vectors ai and ai. The Hamiltonian is written as 

Hbhz = 2j 2j 2j ^^^L.,Q«.^z 

/ a-s,p Sz-± 

-ZZ Z Z'-CUC.V...(86) 

i a,p ii=±x,±y Sz-± 

where a^p are s ox p bands, s^ is the z component of the spin, 
andyu specifies the four nearest-neighbor bonds. The hopping 
term is a matrix in the band basis: 



f;yv_ — 



tspe-''^^^ 



tspe''^^^ 



-t 



pp 



(87) 



where ^^ is the angle between x axis and the direction of the 
bondyu, which can be 0, |, tt, y . 

Fourier transformation of the first term in Eq. (86) gives 

(88) 
Similarly, the hopping term reduces in the k representation to 



-E 



a/3 ^t r^ 



Yj ^k[(^^^ - W) Yj^^^^ k • a^) / / 
k // 

+ (tss + tpp) ^(cos k • a^) o-^ (8) / 
p 

-\- (2tsp sin k • ai) (Ty / 



-\- (Itsp sin k • ai) (Tx <S> s^ 



Ck . 



(89) 
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Therefore, by defining the following Dirac T matrices 



T^ = ..^s, = [ ; S 



T = CTr^S) s. = 



s, 

Sy 

s. 



-// 



r = 0-3; / = 

the total Hamiltonian Eq. (86) is written in the form 



// / ' 
/ 



/f(k) = Jo(k)/ + ^4(k)P 



(90) 
(91) 
(92) 
(93) 
(94) 

(95) 



a=l 



with the coefficients 

Jo(k) = (6, + €p)/2 

-(tss - tsp)(cos k • ai + COS k • ai), (96) 

^i(k) = 0, (97) 

^2(k) = 0, (98) 

J3(k) = 2^,^sink.a2, (99) 

J4(k) = 2^,psink-ai, (100) 

J5(k) = (€s - €p)/2 

-Oss + tsp)(cos k • ai + cos k • ai). (101) 
The energy eigenvalues are given by 



E(k) = do(k)± /^4(k)2. 



(102) 



Note that, since ^-orbital is parity even and /7-orbital is par- 
ity odd, the inversion operator in this band basis is constructed 
as 



n = o-,(8)/ = r\ 



(103) 



One may notice that the T matrices defined above have the 
following symmetry: 



©r©" 



nrn- 



-r^ (a = 1,2,3,4) 

+r« (a = 5) 

-r^ (a = 1,2,3,4) 

+r^ (a = 5) 



(104) 



(105) 



Since only T^ is even under TR and inversion, at a TRIM A/ 
where the system preserves both TRS and IS, the Hamiltonian 
must have the form 



H(k = Ad = do(At)I^ds(At)r 



(106) 



Also, Kramers theorem dictates that l^-, T) and l^-, i) are degen- 
erate at A/, and so are \px + ipy, T) and \px - ipy,i). Remem- 
bering that the s- and /7-orbital bands have even (-\-) and odd 
(-) parities, respectively, we denote those degenerate states at 
Ai as |-h) and |-). One can easily see that the eigenvalues of 
H(Ai) given by Eq. (106) for those states are 

(+|if(A,)|+) = Jo(A/) + ^5(A/) = ^+ 
{-\H(Ai)\-} = do(Ai) - dsiAi) ^ E., (107) 

because <+|r5|+) = 1 and (-ir^h) = -1 due to Eq. (103). 

When dsiAi) > at a TRIM A,, E+ > E. and the |-) state 
is occupied at that A/ if the system is half-filled; in this case, 
the parity of the occupied state at A/, ^(A^) in Eq. (84), is -1. 
On the other hand, when dsiAt) < 0, E+ becomes smaller than 
E- and hence the |-h) state is occupied, making S(Ai) = -hi. 
Therefore, in the present 4-band model one may conclude 



^(A,) = -sgn[J5(A,)]. 



(108) 



If we set |ai| = lail = n, the four TRIMs are (0,0), (0, 1), 
(1,0), and (1, 1). Putting Eq. (101) into Eq. (108), one obtains 



S(A(nun2)) = -sgn 



(109) 
This leads to the following conclusion regarding the topology 
of this 4-band system: 

• If €s - ep > 4(tss -\- tsp), S < for all A/, which leads to 
the Z2 invariant v given by Eq. (85) to become 0. This 
means that the system is a trivial insulator. 



If6v 



< ^(ts 



■ ^p < ^{tss + tsp)^ ^ > at k = but ^ < for other 
three TRIMs, which leads to v = 1. This means that the 
system is a TI. 

It is instructive to note that in the trivial state the band or- 
der is E- < E+, which means /7-orbital bands lie below the 
^--orbital bands, which is normal for a band insulator. When 
the system becomes a TI, the band order flips only at the T 
point in the 2D BZ, and such a band inversion occurring at an 
odd number of TRIMs in the BZ is the key to realizing the 
nontrivial Z2 topology. In the case of HgTe, this band inver- 
sion occurs due to strong SOC. Therefore, this example tells 
us that a viable strategy for discovering a TI is to look for 
materials having strong SOC that causes band inversion at an 
odd number of TRIMs in the BZ. 

4. TI Materials Discovered to Date 

Table I gives a summary of the TI materials discussed in 
this section. This table lists only those materials that have 
been experimentally addressed, since there are too many ma- 
terials that have been predicted to be TIs but not tested exper- 
imentally. 

4. 1 Two-dimensional TIs 

The first material that was experimentally identified as 
a TR-invariant TI was CdTe/HgTe/CdTe quantum well,^^^ 



11 



J. Phys. Soc. Jpn. 



INVITED REVIEW PAPERS 



Table I. Summary of topological insulator materials that have bee experimentally addressed. (In this table, S.S., P.T., and S.M. stand for surface state, phase 
transition, and semimetal, respectively.) 



Type 



Material 



Band gap 



Bulk transport 



Remark 



Reference 



2D, y = 1 


CdTe/HgTe/CdTe 


< 10 meV 


insulating 


high mobility 


26) 


2D, y = 1 


AlSb/InAs/GaSb/AlSb 


~4meV 


weakly insulating 


gap is too small 


64) 


3D(1;111) 


Bii_^Sb;, 


< 30 meV 


weakly insulating 


complex S.S. 


31,35) 


3D(1;111) 


Sb 


semimetal 


metallic 


complex S.S. 


34) 


3D(1;000) 


Bi2Se3 


0.3 eV 


metallic 


simple S.S. 


79) 


3D(1;000) 


Bi2Te3 


0.17 eV 


metallic 


distorted S.S. 


80,81) 


3D(1;000) 


Sb2Te3 


0.3 eV 


metallic 


heavily p-type 


82) 


3D(1;000) 


Bi2Te2Se 


-0.2 eV 


reasonably insulating 


pjcx up to 6 Qcm 


96,99,101) 


3D(1;000) 


(Bi,Sb)2Te3 


< 0.2 eV 


moderately insulating 


mostly thin films 


168) 


3D(1;000) 


Bi2-;cSb;,Te3-3;Se3; 


< 0.3 eV 


reasonably insulating 


Dirac-cone engineering 


103, 104, 187) 


3D(1;000) 


Bi2Tei.6Si.4 


0.2 eV 


metallic 


^-type 


185) 


3D(1;000) 


Bii.iSbo.9Te2S 


0.2 eV 


moderately insulating 


P;cxupto0.1 Dcm 


185) 


3D(1;000) 


Sb2Te2Se 


7 


metallic 


heavily /7-type 


96) 


3D(1;000) 


Bi2(Te,Se)2(Se,S) 


0.3 eV 


semi-metallic 


natural Kawazulite 


186) 


3D(1;000) 


TlBiSe2 


-0.35 eV 


metallic 


simple S.S., largest gap 


87-89) 


3D(1;000) 


TlBiTe2 


-0.2 eV 


metallic 


distorted S.S. 


89) 


3D(1;000) 


TlBi(S,Se)2 


< 0.35 eV 


metallic 


topological P.T. 


93,94) 


3D(1;000) 


PbBi2Te4 


-0.2 eV 


metallic 


S.S. nearly parabolic 


106, 109) 


3D(1;000) 


PbSb2Te4 


7 


metallic 


/?-type 


106) 


3D(1;000) 


GeBi2Te4 


0.18 eV 


metallic 


^-type 


96-98) 


3D(1;000) 


PbBi4Te7 


0.2 


metallic 


heavily ^-type 


110) 


3D(1;000) 


GeBi4Te7 


7 


7 


no data pubHshed yet 


111) 


3D(1;000) 


(PbSe)5(Bi2Se3)6 


0.5 eV 


metallic 


natural hetero structure 


114) 


3D(1;000) 


(Bi2)(Bi2Se2.6So.4) 


semimetal 


metallic 


(Bi2)„(Bi2Se3)m series 


112) 


3D(1;000) 


(Bi2)(Bi2Te3)2 


7 


7 


no data pubHshed yet 


111) 


3D TCI 


SnTe 


0.3 eV (4.2 K) 


metallic 


Mirror TCI, hm = -2 


54) 


3D TCI 


Pbi_;,Sn;,Te 


< 0.3 eV 


metallic 


Mirror TCI, hm = -2 


140) 


3D TCI 


Pbo.77Sno.23Se 


invert with T 


metallic 


Mirror TCI, hm = -2 


138) 


3D(1;111)? 


SmB6 


20meV 


insulating 


possible Kondo TI 


118-121) 


3D (0;001)? 


Bii4Rh3l9 


0.27 eV 


metallic 


possible weak 3D TI 


123) 


3D(1;000)? 


7?BiPt (R = Lu, Dy, Gd) 


zero gap 


metallic 


evidence negative 


130) 


Weyl S.M.? 


Nd2(Iri_;,Rh;,)207 


zero gap 


metallic 


too preliminary 


135) 



namely, a thin layer of HgTe sandwiched by CdTe. This is a 
2D system where the degree of freedom for the perpendicular 
direction is quenched due to the quantum confinement of the 
electronic states in the HgTe unit and the resulting subband 
formation. Both HgTe and CdTe crystallize in zinc-blende 
structure, and their superlattices have been actively studied 
because of their application to infrared detectors. Therefore, 
the necessary technology to synthesize the required quantum 
well was already developed before the discovery of its TI na- 
ture, although it involves a very specialized molecular-beam 
epitaxy (MBE) technique to deal with mercury. ^^^ 

As we have seen in the analysis of the BHZ model in the 
previous section, the inversion between p- and ^--orbital bands 
is essential for the system to obtain the TI nature. Bulk HgTe 
realizes such a band inversion, while CdTe does not. There- 
fore, HgTe is a good starting material for conceiving a TI 
phase; however, there is a problem in the band structure of 



bulk HgTe, that is, a crystal- symmetry-protected degeneracy 
at the r point makes the system to be intrinsically a zero-gap 
semiconductor,^^^ which means that there is no band gap be- 
tween the p- and ^--orbital bands to protect the occupied state 
and makes the topological invariant well defined. However, 
by sandwiching HgTe by CdTe, which has slightly larger lat- 
tice constant, the epitaxial strain exerted on HgTe breaks the 
cubic lattice symmetry and leads to a gap opening, and hence 
the system can become a genuine insulator. It was predicted 
by BHZ that above a certain critical thickness the strained 
HgTe unit retains the band inversion and should be a TI.^^^ 
This prediction was confirmed by a group at the University 
of Wiirzburg lead by Molenkamp via transport experiment of 
micro-fabricated samples. ^^^ They found that, when the thick- 
ness of the HgTe unit is above the critical thickness dc of 6.3 
nm, their samples show a "negative energy gap" (i.e. band is 
inverted) and quantization of the conductance to 2e^/h in zero 
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magnetic field was observed when the chemical potential is 
tuned into the gap. In contrast, when the thickness is below 
dc, the band inversion is lost and they observed diverging re- 
sistivity. 

The Wiirzburg group later reported the observation of non- 
local transport,^^^ which gave further support to the existence 
of an edge state. They also showed that the ID edge state re- 
sponsible for the quantized transport in the 2D TI phase is 
likely to be helically spin polarized, by fabricating an elab- 
orate device structure which relies on the conventional spin 
Hall eff'ect that occurs in metallic (doped) HgTe.^^^ 

The CdTe/HgTe/CdTe quantum well can be made very 
clean with the carrier mobility of up to ~10^ cm^/Vs, which 
makes it possible to study quantum transport properties. On 
the other hand, the drawback of this system is that the bulk 
band gap that opens due to the epitaxial strain is very small 
(up to ~ 10 meV, depending on the thickness^^^), which makes 
the detection of the TI phase to be possible only at very low 
temperature. Also, since the synthesis of CdTe/HgTe/CdTe 
quantum wells requires dedicated MBE machines, the sources 
of samples for basic physics experiments are currently very 
limited, which has made the progress of experimental studies 
of 2D TIs (QSH insulators) relatively slow. 

Recently, another 2D TI system, AlSb/InAs/GaSb/AlSb 
quantum well, was theoretically predicted^^^ and experimen- 
tally confirmed.^"^'^^^ The essential workings of this system 
is the following: The valence-band top of GaSb lies above 
the conduction-band bottom of InAs. Hence, when InAs and 
GaSb are in direct contact and they are both quantum con- 
fined (by the outer units of AlSb which has a large band gap 
and works as a barrier), the resulting hole subband in GaSb 
may lie above the electron subband in InAs, and therefore the 
band order of this quantum well is inverted. The band gap in 
this quantum well arises from anti-crossing of the two sub- 
bands at finite momentum and hence is very small (~4 meV), 
which makes clean observation of the helical edge state very 
difficult. ^"^'^^^ The best evidence for the TI phase has been ob- 
tained via observation of the 2e^ Ih quantization of the zero- 
bias Andreev reflection conductance through Nb point con- 



tacts 



65) 



Other possible candidates of 2D TIs include Bi bilayer,^^^ 
Na2lr03,^^^ graphene with artificially enhanced SOC,^^^ and 
silicene.^^^ Since the experimental studies of 2D TIs have 
been hindered by the scarcity of samples, new discovery of 
2D TI materials are strongly call for. 

4.2 Three-dimensional TIs 

As already mentioned in Sec. 2, the first 3D TI material 
that was experimentally identified was Bii_;cSb;c.^^^ following 
the very specific prediction by Fu and Kane.^^^ This material 
is an alloy between Bi and Sb and it naturally possesses the 
two essential features, (i) band inversion at an odd number of 
TRIMs and (ii) opening of a bulk band gap, in the Sb concen- 
tration range of 0.09 to 0.23.^^^ The 3D Z2 invariant has been 



identified as (1;111). 

Unfortunately, this system turned out to be not very suit- 
able for detailed studies the topological surface state due to 
the complicated surface-state structure,^ ^'^^^ as can be seen 
in Fig. 2. This is because its parent material, Bi, already 
harbors prominent spin-non-degenerate surface states due to 
the strong Rashba eff'ect on the surface of this material.^^^ 
In Bii_;cSb;c. the non-topological, Rashba-split surface states 
are responsible for 2 or 4 Fermi-level crossings of the sur- 
face states (depending on the chemical potential),^ ^'^^^ and 
the topological one contribute just one additional Fermi-level 
crossing. First-principle calculations of the surface states of 
Bii_;cSb;c have been reported,^^'^"^^ but the predicted surface- 
state structure does not really agree with the experimental re- 
sults. Such an incomplete understanding of the nature of the 
surface state is partly responsible for confusions occasionally 
seen in interpretations of experimental data. For example, in 
the STS work which addressed the protection from backscat- 
tering in the surface state in Bii_;cSb;c,^^^ the analysis consid- 
ered the spin polarizations of only those surface states that 
are also present in topologically trivial Bi,^^^ and yet, it was 
concluded that the result gives evidence for topological pro- 
tection. 

Nevertheless, Bii.^^^Sb;^ is unique among known 3D TI ma- 
terials in that it has an intrinsically high 2D carrier mobility of 
-10"^ cm^Vs (despite the fact that it is an alloy), which makes 
it easy to study novel 2D quantum transport. ^^'^^^ Also, it is 
relatively easy for this system to reduce the bulk carrier den- 
sity to -10^^ cm"^ in high-quality single crystals, making it 
possible, for example, to perform Landau level spectroscopy 
of the surface states via magneto-optics. ^^^ The bulk band gap 
of Bii_;cSb;c is not very large (up to ~30 meV depending on 
x)^^^ but it is large enough to detect the 2D transport proper- 
ties at 4 K. 

Along with Bii-j^Sbx, Fu and Kane predicted^^^ that HgTe, 
Qf-phase of Sn (called "gray tin"), and Pbi-jcSn^^Te would be- 
come 3D TIs under uniaxial strain to break the cubic lattice 
symmetry. They also suggested that Bi2Te3 would be a candi- 
date, but they did not perform band calculations to elucidate 
the parity eigenvalues. Such calculations were done by Zhang 
et al,^^'' who came up with a concrete prediction that Bi2Se3, 
Bi2Te3, and Sb2Te3 should be 3D TIs but Sb2Se3 is not; fur- 
thermore, Zhang et al proposed a low-energy eff'ective model 
to describe the bulk band structure of this class of materials. 
This model, with some corrections made later,^^^ has become 
a popular model for theoretically discussing the properties of 
a 3D TI. Experimentally, existence of a single Dirac-cone sur- 
face state was reported in 2009 for Bi2Se3 by Xia et aC^'' and 
for Bi2Te3 by Chen et al}^^ and also by Hsieh et al ;^^^ Sb2Te3 
was measured by Hsieh et al}^'' along with Bi2Te3, but the ex- 
istence of the topological surface state was left unconfirmed 
due to the heavily /7-type nature of the measured samples. The 
topological nature of Sb2Te3 was confirmed only recently in 
thin-film samples using STS.^^^ 
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Fig. 6. (Color online) Crystal structure of tetradymite chalcogenides, 
Bi2Se3 (Xl =X2 = Se), Bi2Te3 (Xl = X2 = Te), and Bi2Te2Se (XI = Te; 
X2 = Se). The quintuple layer (QL) enclosed by brown cage is the build- 
ing block of this type of materials. The stack of the QLs are in the -A-B-C- 
A-B-C- manner, and hence the unit cell consists of a stack of 3 QLs. The 
Bi2-;cSb;cTe3_3;Se3; (y > 1) compound also takes this structure, in which X2 = 
Se, XI = Se/Te, and the original Bi site becomes a mixture of Bi/Sb. Taken 
from Ref. 99; copyright American Physical Society (2010). 



Bi2Se3, Bi2Te3, and Sb2Te3 all crystallize in tetradymite 
structure, which consists of covalently bonded quintuple lay- 
ers (e.g. Se-Bi-Se-Bi-Se) that are stacked in -A-B-C-A-B-C- 
manner and are weakly interacting with van der Waals force 
(Fig. 6); therefore, those materials cleave easily between quin- 
tuple layers (QLs). Since each QL is about 1 nm thick, the 
lattice constant along the c-axis is about 3 nm. 

The 3D Z2 invariant of these tetradymite systems is (1;000), 
which means that topological Dirac-cone surface state is cen- 
tered at the r point of the surface BZ (Fig. 7). This simplic- 
ity of the topological surface state and the absence of non- 
topological surface states make those materials well suited for 
experimentally addressing the properties of the topological 
surface state. Also, single-crystal growth of those materials is 
relatively simple, which made it easy for many experimental- 
ists to start working on them. Furthermore, the bulk band gap 
of Bi2Se3 is relatively large, 0.3 eV, and thus one can see tech- 
nological relevance that topological properties of this material 
may potentially be exploited at room temperature. All those 
factors helped initiate a surge of research activities on 3D TIs. 

The surface-state structure of Bi2Se3 is relatively simple 
and presents an almost idealized Dirac cone with only slight 
curvature, as shown schematically in Fig. 7(a). In contrast, 
the surface state of Bi2Te3 is a bit more complicated [see Fig. 
7(b)] and the Dirac point is located beneath the top of the va- 
lence band, which makes it difficult to probe the surface trans- 
port properties near the Dirac point without being disturbed 
by bulk carriers in Bi2Te3. 

Another difference between the two materials is that the 
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Fig. 7. (Color online) Schematic bulk and surface band structures of (a) 
Bi2Se3 and (b) Bi2Te3. Note that the surface states are spin non-degenerate 
and are helically spin polarized. Representative constant-energy contours of 
the Dirac cones for (c) Bi2Se3 and (d) Bi2Te3 are also schematically shown. 
Note that the spin vector is always perpendicular to the wave vector ki in both 
Bi2Se3 and Bi2Te3, but the Fermi velocity vector \f can be non-orthogonal 
to the spin vector in Bi2Te3 due to the hexagonal warping, which leads to 
strong quasiparticle interference. 



constant-energy contour of the Dirac cone is almost spherical 
in Bi2Se3 [Fig. 7(c)], while it presents significant hexagonal 
warping in Bi2Te3 [Fig. 7(d)]. This warping is caused by a P 
term which stems from cubic Dresselhaus spin-orbit coupling 
at the surface of rhombohedral crystal sy stems. ^^^ Intrigu- 
ingly, this hexagonal warping gives rise to peculiar physics 
such as strong quasiparticle interference^^' ^^^ and appearance 
of a finite out-of-plane spin polarization. ^^^ 

The discovery of tetradymite TI materials demonstrated 
that elucidating the parity eigenvalues based on ab initio band 
calculations is a practical and powerful strategy for making 
predictions for new TI materials. In fact, the next TI ma- 
terial discovered after Bi2Se3 and Bi2Te3 was TlBiSe2, and 
its discovery, ^^~^^^ done in 2010, was also lead by theoretical 
predictions.^^' ^^^ By 2010, the competition in the TI research 
was already very heated and two theoretical groups indepen- 
dently made predictions that thallium-based III-V-IV2 ternary 
chalcogenides T1M'X2 (M' = Bi and Sb; X = S, Se, and Te) 
should be TIs; experimentally, the TI nature of TlBiSe2 was 
reported first by Sato et al.,^^'' which was soon followed by 
independent reports by Kuroda et al}^^ and by Chen et al}"^^ 
The latter also confirmed TlBiTe2 to be a TI.^^^ The surface 
state structure of TlBiSe2 is similar to that in Bi2Se3, and its 
simplicity makes it suitable for studying fundamental proper- 
ties of the Dirac cone without being bothered by additional 
complications. Also, Sato et al elucidated that the bulk band 
gap of TlBiSe2 is 0.35 eV,^^^ which is among the largest so 
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far reported for 3D TI materials. 

Intriguingly, TlBiS2, which was initially predicted to be a 
jj 90,91) turned out to be non-topological. This indicates that 
ab initio calculations are not always accurate in every de- 
tails, and one should interpret the results with due care.^^^ 
The non-topological nature of TlBiS2 naturally points to a 
topological phase transition occurring in the solid solution 
TlBi(Si_;cSe;c)2- Such a transition was studied independently 
by Xu et al^^^ and by Sato et al^^'' and was found to be lo- 
cated at X = 0.5. Intriguingly, Sato et al found an unexpected 
gap opening at the Dirac point when the composition is on the 
topological side and is close to the transition.^"^^ This is sur- 
prising, because such a gap opening points to a lifting of the 
Kramers degeneracy, but in this system TRS is not explicitly 
broken. While the origin of this "Dirac gap" is not clear and 
its existence has been a matter of debate, a recent follow-up 
experiment substantiated its unique properties. ^^^ 

After the discovery of TlBiSe2, various ternary compounds 
have been identified to be 3D TIs. In a preprint that was never 
published, Xu et al reported ARPES data to show topolog- 
ical surface states in GeBi2Te4, Bi2Te2Se, and Sb2Te2Se.^^^ 
In this preprint, Bi2Te2Se and Sb2Te2Se were reported to be 
n- and /?-type, respectively, while GeBi2Te4 was suggested 
to be an intrinsic insulator; this result on GeBi2Te4 was pub- 
lished later in a diff'erent paper. ^^^ However, a follow-up work 
by Okamoto et al showed that GeBi2Te4 is naturally electron 
doped.^^^ The GeBi2Te4 crystals grown in the present author's 
laboratory were also highly conducting with ^-type carriers. 

The story about Bi2Te2Se is interesting. While Xu et 
al reported metallic ^-type nature for the stoichiomet- 
ric compound,^^^ Ren et al discovered^^^ that by grow- 
ing crystals from a slightly Se-rich starting composition of 
Bi2Te1.95Se1.05, one can obtain crystals showing a large bulk 
resistivity exceeding 1 Dcm. Furthermore, they demonstrated 
that in such crystals, the chemical potential is located within 
the bulk band gap and one can observe clear SdH oscillations 
coming from the topological surface state, which contributed 
~6% of the total conductance in a 260-yum-thick bulk crys- 
tal. ^^^ This value is to be contrasted with the preceding surface 
transport study by Qu et al on Bi2Te3, which found the sur- 
face contribution of only ~0.3% in a 100-yum-thick sample. ^^^^ 
In fact, Bi2Te2Se was the first 3D TI material to present a rea- 
sonably bulk-insulating behavior, which opened the door for 
detailed transport studies of the topological surface state. This 
discovery by Ren et al was followed by an independent report 
by Xiong et al}^^^ who reported an even larger resistivity of 
6 Dcm. Detailed defect chemistry in Bi2Te2Se was discussed 
hy YmetaO^^^ 

Ren et al tried to further improve the bulk-insulating prop- 
erty of Bi2Te2Se, and they reported two possible routes. 
One is to expand the phase space of the compositions into 
Bi2-;cSb;cTe3_3;Se3;, and they found a series of particular com- 
binations of (x, y) where the samples show maximally bulk- 
insulating behavior. (This material is discussed in detail in 



Sec. 7.) The second route to improve the property of Bi2Te2Se 
is to employ Sn doping to the Bi site.^^^^ Using both routes, it 
has been shown to be possible to achieve surface-dominated 
transport in bulk single crystals. ^^^' ^^^^ 

Another relatively new TI compound is Pb(Bii_;cSb;c)2Te4, 
for which the observation of a topological surface state, as 
well as a switching between p- and ^-type surface carriers 
with X, has been reported by Souma et al}^^^ which was mo- 
tivated by theoretical predictions.^^' ^^^' ^^^^ Kuroda et al in- 
dependently reported the topological surface state for one of 
its end members, PbBi2Te4.^^^^ Interestingly, Pb-based mate- 
rials form various homologous series of compounds, such as 
^PbTe-mBi2Te3. In fact, PbBi2Te4 can be viewed as a mem- 
ber of this homologous series with n = m = 1 . Observation of a 
topological surface state was also reported for its ^ = 1 , m =2 
member, PbBi4Te7.^^^^ In passing, Ge-based materials form 
similar homologous series as Pb-based ones. GeBi2Te4^^'^^^ 
and GeBi4Te7^^^^ are examples with ^ = m =1 and ^ = 1, m 
=2, respectively. 

It is known that Bi bilayer and Bi2Se3 (or Bi2Te3) also 
form a homologous series of compounds, (Bi2)„(Bi2X3)^ 
iX = Se, Te).^^^^ Topological surface state has been re- 
ported for a partially-S-substitued variant, (Bi2)(Bi2Se2.6So.4) 
(= Bi4Se2.6So.4) which is a semimetal.^^^^ In this type of ho- 
mologous series, (Bi2)(Bi2Te3)2 (= BiTe) has also been syn- 
thesized and was found to be topological^ ^^^ 

As mentioned in conjunction with Bii_;cSb;c, one of the ini- 
tial predictions of Fu and Kane was that strained HgTe would 
be a 3D TI with the 3D Z2 invariant (1;000).^^^ This was 
experimentally verified in 2011 using MBE-grown films by 
Brtine et aO^^^ who reported that a 70-nm-thick HgTe film 
grown on a CdTe substrate is uniformly strained (0.3%) and 
obtains a band gap of ~22 meV; a strange pattern of quantum 
Hall eff'ect was observed with Hall plateaus appearing at fill- 
ing factor y = 2, 3, 4, 5, 7, 9. This was interpreted^ ^^^ to be a 
result of a combination of two quantum Hall eff'ects occurring 
in parallel at the top and bottom surfaces that contain slightly 
diff'erent carrier densities. So far, strained HgTe is unique in 
that the mobility of the surface carriers can be high enough to 
present quantum Hall eff'ect. 

An interesting material which was experimentally found 
to be topological without any theoretical prediction is the 
(PbSe)5(Bi2Se3)3^ compound (m = 1, 2).^^"^^ This material 
forms a natural multilayer hetero structure consisting of topo- 
logical insulator (Bi2Se3) and an ordinary insulator (PbSe) 
units (Fig. 8). The precise atomic positions within the unit cell 
has not been determined yet for this material, and hence no 
ab initio band calculations have been performed. Nakayama 
et al found that in this material the PbSe unit provides a bar- 
rier for the electronic states in the Bi2Se3 unit to be quantum 
confined, which, for m = 2, leads to hybridization of the topo- 
logical interface states at the top and bottom of each Bi2Se3 
unit, resulting in a gapped Dirac-cone structure. Also, the bulk 
band gap of this m = 2 compound was found to be as large as 
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Fig. 8. (Color online) (a) Schematic picture of the crystal structure of 
(PbSe)5(Bi2Se3)3m, which forms a natural multilayer heterostructure consist- 
ing of topological insulator (Bi2Se3) and an ordinary insulator (PbSe) units, 
(b) ARPES intensities of (PbSe)5(Bi2Se3)3m for m = 1 (left) and 2 (center), 
together with that of Bi2Se3 (right), plotted as a function of the binding en- 
ergy and wave vector measured along the ky axis (y - K cut) at T = 30 K 
with hy = 60 eV. (c) Magnified view of the same data for m = 1 (left) and 2 
(center) at low energy; the right panel shows the second derivative of the m 
= 2 data, in which the gap opening at the Dirac point is evident. Taken from 
Ref. 1 14; copyright American Physical Society (2012). 



0.5 eV, which is probably due to the quantum confinement ef- 
fect on the bulk band. For m = 1 , the quantum confinement 
apparently lifts the non-trivial Z2 topology from the Bi2Se3 
unit and gives rise to a degenerate parabolic band (Fig. 8). 

4.3 Candidate 3D TIs 

There are other candidate materials of 3D TIs for which the 
topological surface states have not been observed. There have 
been many predictions for 3D TIs based on ab initio band 
calculations/ ^^^ but in the following I will discuss only those 
materials that have been experimentally addressed. 

A very interesting candidate is SmB6, which may be a topo- 



logical Kondo insulator, ^^^'^^^"^ because the band bap of this 
material stems from Kondo eff'ect. If this material is indeed a 
TI, it will be the first material where electron correlation and 
nontrivial band topology both play important roles. So far, 
strong evidence for surface-dominated transport at low tem- 
perature has been reported/ ^^~^^^^ but existing ARPES data 
do not resolve topological surface state within the small gap 
(~20 meV).^^^^ A puzzling feature in the reported surface- 
transport data^^^^ is that the observed surface carrier density 
of 1.1 X 10^^ cm"^ is about 1000 times larger than that found 
in other TI materials and is obviously too large for the ex- 
pected Dirac carrier density. ^^^^ One should keep in mind that 
the surface-dominated transport can be observed in insulators 
in ambient condition due to the formation of a trivial accumu- 
lation (or inversion) layer, as reported for pure Te by von Kl- 
itzing and Landwehr,^^^^ and hence such an observation alone 
does not give evidence for a TI phase. 

Another interesting candidate is Bii4Rh3l9, which may be 
a weak 3D TI with the Z2 index (0;001).^^^^ The crystal struc- 
ture of this material can be viewed as a stack of intermetallic 
[(RhBi4)3l]^^ layers and intervening [Bi2l8]^~ zigzag-chain 
layers; the former consist of graphene-like honeycomb net- 
work of edge- sharing RhBig cubes, giving rise to an electronic 
structure which is similar to that in graphene but is enriched 
by strong SOC due to Bi. The ab initio band calculations pre- 
dict that this is a quasi-2D system which essentially consists 
of a stacked 2D TI layers; as a result, the system is predicted 
to be a weak 3D TI with the topological surface state expected 
only on the sides of the stacks. The bulk band structure seen 
by ARPES is in reasonable agreement with the calculations, 
but no surface state was detected because the ARPES was 
done on the cleaved top surface where the surface state is not 
expected. Obviously, some transport measurements to detect 
the topological surface state on the side surface is necessary 
to elucidate the TI nature of this material. 

4.4 Topological semimetal 

The term "topological semimetal" can be used for three dif- 
ferent classes of materials. The first is an ordinary semimetal 
possessing separate electron and hole pockets compensating 
each other, but the parity eigenvalues of the valence band 
gives rise to a nontrivial Z2 topology; an example is pure 
Sb, whose topology is characterized by the 3D Z2 invari- 
ant (l;!!!).^^'^"^^ The second is a zero-gap semiconductor 
which possess band inversion, such as unstrained HgTe; this 
class of materials can be made genuine TIs by lowering the 
crystal symmetry to open a gap. The third is called Weyl 
semimetal,^^^^ which is probably the most fundamentally in- 
teresting. In the following, I will discuss the latter two classes. 

Interesting examples of the candidate topological zero-gap 
semimetals are Heusler or half-Heusler compounds. ^^^^ They 
are ternary intermetallic compounds with X2YZ (Heusler) 
or XYZ (half-Heusler) compositions. Their crystal structures 
consist of basic YZ sublattice (which takes the zinc-blende 
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Structure similar to HgTe) stuffed with X in the void space. 
Many Heusler or half-Heusler compounds were predicted to 
have band inversion and be topological based on ab initio 
band calculations ;^^^~^^^^ under uniaxial pressure which low- 
ers the crystal symmetry, they obtain a gap and are expected to 
become TIs. The ARPES experiments done on 7?PtBi (R = Lu, 
Dy, Gd) found metallic surface states distinct from the bulk 
bands, but the surface-band structure does not appear to sup- 
port a topological origin. ^^^^ Interestingly, superconductivity 
has been found in LaPtBi^^^^ and YPtBi;^^^^ if those materi- 
als are indeed topological, the fate of the surface states in the 
superconducting state is an intriguing issue. 

The idea of Weyl semimetal was theoretically proposed by 
Wan et al. in 2011.^^^^ Weyl fermions have definite chiral- 
ity, which protects them from gapping, and they are described 
by a massless two-component Dirac equation. In solids, when 
TRS is broken in an inversion symmetric system and non- 
degenerate valence and conduction bands touch at an acciden- 
tal degeneracy point in a 3D BZ, the low-energy physics is ap- 
proximated by massless two-component Dirac equation and 
the quasiparticles obey the Weyl Hamiltonian H = ±vfO- • k. 
By symmetry, such a degeneracy point (Weyl point) must 
come in pairs, and Wan et al showed that a pair of Weyl 
points give rise to an arc of zero-energy excitation (Fermi arc) 
to connect them in the projected surface BZ.^^^^ The appear- 
ance of this gapless state on the surface of Weyl semimetal 
can be viewed as a result of the bulk-edge correspondence 
to signify the non-trivial topology. ^^^'^^"^^ It was proposed 
based on ab initio band calculations that pyrochlore iridates 
such as Y2lr207 in the antiferromagnetic phase may realize 
such a Weyl semimetal. ^^^^ In this context, a recent study of 
Nd2(Iri_;cRh;c)207 found a transition from a correlated metal 
to a Mott insulator with decreasing x that is accompanied by 
a gradual reduction of the spectral weight, and it is possible 
that the Weyl semimetal phase exists at the critical point. ^^^^ 

The Weyl semimetal phase is also conceivable in a TR sym- 
metric system with broken inversion symmetry, and it was 
proposed that such a phase may be achieved at the topo- 
logical phase transition point of TlBi(Si_;cSe;c)2 by breaking 
the inversion symmetry through an elaborate multilayer struc- 
ture. ^^^^ 

4.5 Topological crystalline insulator 

We have already seen that band insulators can be topolog- 
ically classified by evaluating the Z2 invariant from valence- 
band Bloch wave functions. This classification is based on 
TRS of the system. It turned out that this is not the only possi- 
ble topological classification of band insulators, but it is also 
possible to classify band insulators based on topologies pro- 
tected by point-group symmetries of the crystal lattice. ^^^^ 
Those insulators that have nontrivial topology protected by 
point-group symmetries are called topological crystalline in- 
sulators (TCIs).^^^^ So far, concrete topological invariants are 
elucidated for systems possessing four-fold (C4) or six-fold 



(Ce) rotation symmetry ^^^^ and also for systems with mirror 
symmetry.^^^ In particular, the latter case gained significant at- 
tention after the prediction of a concrete example, SnTe, was 
made by Hsieh et alP'^ 

In TCIs with mirror symmetry, the topology is specified 
by the topological invariant nM called mirror Chern number, 
which evaluates the Chern number in only one of the two 
Hilbert subspaces divided according to the mirror eigenval- 
ues. ^^^ Based on ab initio calculations, Hsieh et al. showed 
that SnTe has band inversions at four TRIMs (four equiva- 
lent L points) in the 3D BZ,^^^ which gives rise to a nontrivial 
mirror Chern number nj^ = -2, while its Z2 invariant is triv- 
ial, (0;000). They further predicted that the surface state on 
the {001} surface should present a peculiar double-Dirac-cone 
structure near the X points, because two L points are projected 
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Fig. 9. (Color online) ARPES data on cleaved (001) surface of SnTe. (a) 
The bulk BZ and the corresponding (001) surface BZ of SnTe. The (110) 
mirror plane is indicated by the shaded area. Note that two L points are pro- 
jected to the same X point, (b) ARPES intensity mapping in the surface BZ 
at Ef measured with hv = 21.2qY at 30 K. (c) NQar-Ef ARPES intensity as 
a function of the wave vector and the binding energy measured along the cut 
which is nearly parallel to the X - M direction and is crossing the A2 point 
[one of the red arrows in panel (b)]. (d) Near-£^i7 ARPES intensity along the 
f -X cut [yellow arrow in panel (b)]. (e) Schematic 2D band dispersions near 
Ef concluded from the data in panels (c) and (d) depicting the characteristic 
double Dirac-cone structure. Taken from Ref. 54. 
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onto the same X point [Fig. 9(a)]. ^^^ 

It is interesting that the experimental verification of this 
TCI prediction showcases how severe the competitions in the 
TI research are. The prediction^^^ was posted to the preprint 
server in February 2012. In June, three preprints reporting 
experimental discovery of a TCI were posted to the preprint 
server in about two weeks, and two of them were published on 
the same day in September 2012 in diff'erent journals (Nature 
Physics^^^ and Nature Materials ^^^^), while the other one^^^^ 
reported by a group at Princeton University was never pub- 
lished; the Princeton group performed new measurements, 
and a paper containing new data were submitted in August 
and published in November. ^^^^ 

Among the first two published papers, the one by Tanaka et 
al^"^^ reported straightforward confirmation of the prediction 
in SnTe, and the observed surface state with the double Dirac- 
cone structure [Fig. 9] was in good qualitative agreement with 
the theory; they also showed that the cousin material PbTe 
does not present any surface state. The other one by Dziawa et 
al}^"^^ reported the TCI phase in Pbo.77Sno.23Se, which shows 
a transition to a trivial phase upon increasing the tempera- 
ture. The other one by Xu et al published later^"^^^ reported a 
TCI phase in Pbo.6Sno.4Te and a trivial phase in Pbo.8Sno.2Te, 
together with spin-resolved ARPES data showing helical po- 
larization on each of the double Dirac cones. A more recent 
paper by Tanaka et al}^^^ nailed down that the topological 
phase transition in Pbi-j^Sn^^Te occurs at Xc - 0.25; further- 
more, they found that the separation between the two Dirac 
cones near the X points systematically narrows when x is re- 
duced toward x^, but they never merge before the transition 
eliminates them. 

In passing, the mirror Chern number riM can also be used 
for TR-invariant 3D TIs to further classify them.^^^ For exam- 
ple, Bii_;cSb;c is a TI with Z2 invariant (1 ;1 1 1), and it can have 
hm = ± 1 . The sign of um is called mirror chirality, which is 
related to the sign of the g factor. The first experimental work 
that addressed this additional topological property in a TI was 
the spin-resolved ARPES done by Nishide et al,^^^ who elu- 
cidated that the mirror chirality is -1 in Bii_;cSb;c- 

The discovery of TCIs significantly widened the scope of 
topological materials. Already, detailed topological classifi- 
cation schemes for all point-group symmetries have been pro- 
posed, ^"^^^ and also the mirror topology has been expanded to 
superconductors. ^"^^"^"^^^ Experimentally, elucidating the inter- 
play between Z2 topology and mirror topology in materials 
like SnTe under uniaxial strain would be an interesting issue. 

5. How to Confirm TI Materials 

In this section, I briefly summarize the possible experimen- 
tal procedure to confirm whether a material is a TI or not. In 
the case of 2D TIs, one needs to probe the existence of helical 
ID edge state, which is possible only through quantum trans- 
port experiments using nano-fabricated device structures. The 
existence of the edge state can be seen through conductance 
quantization in the insulating regime.^^^ Also, the helical spin 



polarization of the edge state may be detected by transport 
experiments using spin Hall eff'ect.^^^ 

For 3D TIs, the simplest and the most convincing is to 
observe the Dirac cone by ARPES experiments. To firm up 
the identification of a TI, one should employ spin-resolved 
ARPES to confirm that the Dirac cone is non-degenerate and 
is helically spin polarized.^^'^^^ 

Unfortunately, not all materials are suitable for ARPES, 
which requires clean and flat surface that is usually obtained 
by cleaving single crystals. When single crystals are not avail- 
able or the material does not cleave well, APRES becomes 
difficult. In such a case, one may rely on transport experi- 
ments. Ideally, if the bulk is sufficiently insulating and the 
surface carriers have high enough mobility, one would be able 
to confirm that the transport is occurring through the surface 
and that the surface carriers are Dirac fermions. The former 
can be done by looking at the sample-size dependence of the 
conductance, ^^"^'^^^"^^^^ and the latter may be accomplished 
by elucidating the n Berry phase in the quantum oscillations 
from the surface state.^^' ^^^' ^^^' ^^^~^^^^ (Detailed discussions 
on the identification of the n Berry phase will be given in 
Sec. 8.3.) It should be emphasized that confirming the Dirac- 
fermion nature of the surface carriers is important, because 
trivial accumulation layer or inversion layer that may form 
on the surface of an insulator may also give rise to surface- 
dominated transport. ^^^^ 

The Dirac-fermion nature may also be confirmed by 
STS experiments in magnetic fields, because massless Dirac 
fermions present peculiar Landau quantization in which the 
level spacing changes as Va^ and the zero-energy Lan- 
dau level is pinned to the Dirac point (see detailed discus- 
sions in Sec. 8.2); by looking at the bias- voltage dependence 
of the Landau quantization peaks, one can identify Dirac 
fermions.^^'^^^'^^^^ Similarly, magneto-optics experiments to 
detect the Landau level transitions can be used for detecting 
the peculiar quantization scheme to identify Dirac fermions 
on the surface. '^^^ 

One should note that in reality, it is often very difficult 
to obtain sufficiently bulk-insulating samples of a candidate 
material. In that case, transport measurements are not very 
useful. If it is possible to detect the equilibrium spin current 
(which is carried only by the helical Dirac fermions and hence 
is not bothered by bulk carriers) by some electromagnetic 
means, ^^^^ it would become a very useful tool for identifying 
a TI. However, feasibility of such an experiment is not clear 
at the moment. 

6. Syntheses of TI Materials 

6. 1 Bulk single crystals 

Except for Bii-^^Sbjc, all the confirmed 3D TI materials are 
chalcogenides (i.e. compounds containing chalcogen atoms S, 
Se, and Te). Since chalcogen atoms are volatile, the syntheses 
of chalcogenides are done in sealed evacuated quartz-glass 
tubes, which can sustain temperatures up to 1000°C. Such a 
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necessity of containment limits the range of applicable growth 
techniques, and one usually uses the Bridgman method. In 
this method, the temperature of the melt is gradually reduced 
while keeping a temperature gradient in the tube, so that the 
crystallization starts at the cold end and the crystal grows as 
the solidification proceeds from this end. Most of the popu- 
lar TI materials including Bi2Se3, Bi2Te3, and Bi2Te2Se are 
grown by the Bridgman method. 

Another possible crystal growth technique is the vapor 
transport. In this technique, one puts a chunk of polycrys- 
talline material on one end of a sealed quartz-glass tube. The 
tube is kept for a long time in a furnace with a certain tem- 
perature gradient, in which the polycrystalline material is on 
the hotter end. The temperature gradient is chosen so that the 
material sublimates at the hotter end and crystalizes at some- 
where in a colder part. When some reagent such as I2 is used 
as a transporter, the technique is called chemical vapor trans- 
port; if the material easily sublimates, one does not need a 
transporter and the technique is called physical vapor trans- 
port. Single crystals of TCI materials, SnTe, (Pb,Sn)Se, and 
(Pb,Sn)Te are usually grown with a vapor transport technique. 

Since the tetradymite compounds cleave easily, they can be 
made into a very thin (down to only a few nm thick) flakes 
by employing an exfoliation technique similar to that used 
for making graphene samples using Scotch tapes. Such thin 
flakes are particularly useful for experiments involving gat- 
ing to electrostatically control the surface chemical poten- 

^j^j 153, 154) 

6.2 Thin films 

For the growth of high-quality epitaxial thin films of TI ma- 
terials, molecular beam epitaxy (MBE) technique is usually 
employed. ^^^^ So far, reports of MBE growths have been made 
for Bii_;cSb;c^^^^ or relatively simple tetradymite compounds 
Bi2Se3,^'^^'^^^'^^^-^^^^ Bi2Te3,i6^> and Sb2Te3,^^^ as wefl as 
their solid-solutions such as (Bi,Sb)2Te3.^^^^ Those materi- 
als can be grown by co-evaporating the constituent elements 
with suitable flux ratios. Chemical vapor deposition (CVD) 
technique has also been applied to the growth of Bi2Se3,^^^^ 
but the reported film quality has not been as good as the best 
MBE-grown films. 

For epitaxial growths of thin films, the lattice matching be- 
tween the substrate and the grown material is usually very 
important. However, in the case of tetradymite TI materials, 
thanks to the existence of the van der Waals gap between the 
QLs, the lattice matching with the substrate is not crucial and 
the epitaxial growth proceeds in the so-called van der Waals 
epitaxy mode,^^^^ in which the substrate and the films are 
only weakly bonded with the van der Waals force and hence 
the lattice matching condition is greatly relaxed. In particular, 
epitaxial growths of Bi2Se3 have been reported for various 
substrates including Si(lll),^^^~^^^^ graphene-terminated 6H- 
SiC(0001),^50'^6^^ GaAs(lll),^62^ sapphire(0001),^49,i63,i64) 
CdS(0001),i65> andInP(lll).i66) 



It turned out that the control of the substrate temperature is 
the most important parameter for obtaining high-quality films 
of Bi2Se3 with large area of atomically flat terraces. ^^"^^ In par- 
ticular, a two-step method in which the depositions of the first 
layer and subsequent layers are done at diff'erent temperatures, 
to promote initial adhesion and crystallization separately, has 
been reported to yield best quality films.^49, 158, 159, 163, 164) 

For the TCI material SnTe, a technique called hot- wall epi- 
taxy has been used in the past,^^^^ yielding good quality sam- 
ples with a reasonably high mobility (-2700 cm^/Vs). For 
SnTe which is a cubic material with rock-salt structure, good 
lattice matching is crucial for epitaxial growth. BaF2 has been 
traditionally used as a substrate, ^^^^ but recently, using Bi2Te3 
as a buff'er layer between SnTe and sapphire was reported 
to yield high quality films that present surface SdH oscilla- 



tions 



172) 



6.3 Nanoribbons and nanoplates 

Mesoscopic transport experiments of Bi2Se3,^^^^ 
Bi2Te3,^^^'^^^^ and Bi2Te2Se^^^^ have been performed 
using nanoribbons and nanoplates. Bi2Se3 nanoribbons are 
usually synthesized by gold-catalyzed vapor liquid solid 
(VLS) technique. ^^^^ Typically, Bi2Se3 powder is placed in 
the center of a tube furnace through which Ar gas flows and 
transports evaporated Bi2Se3; a silicon substrate decorated 
with Au nanoparticles are placed downstream in the tube 
furnace, and the growth of nanoribbons proceeds as Au 
nanoparticles absorb Bi2Se3 vapor and leave crystallized 
Bi2Se3 nanoribbons beneath them. Naturally, the diameter 
of grown nanoribbons is roughly determined by the size of 
the Au nanoparticles (typically ~20-nm diameter). When 
there are no Au nanoparticles to work as catalyst, nanoplates 
(typically a few nm thick and a few yum wide) are obtained 
instead of nanoribbons. ^^^' ^^^^ 

7. Defect Chemistry and Engineering 

An important theme in the research of TI materials is to 
reduce unintentionally-doped bulk carriers that hinder obser- 
vations of surface transport properties. Only in the case of 
HgTe thin films, which can be grown in very high purity us- 
ing an MBE technique, unwanted bulk carriers are not an is- 
sue. In most other materials, one needs to find suitable ways 
to reduce bulk carriers. In this section, I will use the Kroger- 
Vink notation to describe defects in crystals; for example, V^* 
means a selenium vacancy with double positive charge, and 
Bi!j,g means a bismuth ion sitting on the tellurium lattice site 
with a single negative charge. 

Naturally-grown Bi2Se3 crystals are always electron doped, 
with the typical bulk carrier density ^3d of 10^^ cm"^,^^^~^^^^ 
because of the thermodynamically inevitable Se vacancies 
(y**) or Se anti-site defects (Se*.) that have low formation 



energies 



180,181) 



Doping Ca^^ to the Bi^^ site has been re- 
ported to be useful for reducing ^-type carriers and eventu- 
ally achieving n-to-p type transition. ^^^^ However, Ca doping 
obviously introduces strong scattering centers and the elec- 
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tron mobility becomes low in Ca-doped Bi2Se3.^^^^ Optimiza- 
tion of the growth condition^^^^ and isovalent Sb doping to the 
Bi site^^^' ^^^^ are both reported to be useful for reducing ^3d 
down to 10^^ cm"^ while keeping a high electron mobility that 
allowed observation of surface SdH oscillations/"^^^ although 
n-to-p type transition was not achieved with these routes. Re- 
cently, Ren et al}^^^ succeeded in growing high-mobility p- 
type crystals of Bi2Se3 by combining Cd doping and a Se-rich 
growth condition; furthermore, subsequent annealing of such 
/7-type crystals to compensate for the /7-type carriers by the 
electrons coming from Se vacancies made it possible to ob- 
tain both n- and /7-type samples with low bulk carrier density, 
in which surface SdH oscillations were observable. ^^^^ 

In contrast to Bi2Se3, Bi2Te3 can be grown in both n- and 
p-types}^^^ This is because Te anti-site defects (Te*.) or Bi 
anti-site defects (Bi^j,^) are formed in Te-rich or Bi-rich condi- 
tions, respectively, and Te* . is a donor while Biij,^ is an accep- 
^Qj, 180, 181) Nevertheless, it is difficult to obtain Bi2Te3 samples 
with low ^3D. By growing a boule of crystal with a composi- 
tional gradient, Qu et ah were able to pick up samples with 
^3D ~ 10^^ cm"^ in which they observed surface SdH oscilla- 
tions.ioo) 

The Bi2Te2Se compound can be much more insulating than 
Bi2Se3 or Bi2Te3 if grown in a slightly Se-rich condition.^^^ 
This is because Bi2Te2Se naturally crystallize in a chalcogen- 
ordered structure with Te-Bi-Se-Bi-Te QL units (Fig. 6).^^"^^ 
This peculiar structure solves the problems of Se- vacancies in 
Bi2Se3 and Bi/Te antisite defects in Bi2Te3 at the same time: 
First, because Se is concealed in the middle of the QL, Se 
vacancy is much more difficult to occur than in Bi2Se3. Sec- 
ond, since the electronegativity of Se is stronger than that of 
Te, Bi is more strongly bonded to Se than to Te, which dis- 
courages the occurrence of Bi/Te antisite defects. As a result, 
Bi2Te2Se can be reasonably bulk-insulating. In addition, its 
ordered structure ensures a high mobility, and consequently 
the topological surface state can be easily studied through 
transport properties of this material. ^^' ^^^^ In passing, materi- 
als with similar chalcogen-ordered structure, Bi2Tei.6Si.4,^^^^ 
Bii.iSbo.9Te2S,^^^^ andBi2(Te,Se)2(Se,S) (natural Kawazulite 
mineral) ^^^^ have been recently reported, and the latter two are 
found to show moderately bulk-insulating behavior. 

Importantly, even higher bulk-insulating properties than in 
Bi2Te2Se have been achieved by utilizing a solid solution 
Bi2-xSb;cTe3-3;Se3; (y > I).i03,i04) j^^ ^j^^^ material, the middle 
of the QL is preferentially occupied by Se, hence retaining 
the essential virtue of Bi2Te2Se. In addition, the Bi/Sb ratio 
in the cation layers and the Te/Se ratio in the outer layers af- 
fect the levels of acceptors and donors, respectively, making 
it possible to achieve a maximally compensated situation. In- 
deed, it was found that a reasonably bulk-insulating behavior 
is observed for a series of compositions (x, y) which forms 
a curved line in the x vs j phase diagram (Fig. 10). Further- 
more, the surface Dirac-cone structure of Bi2-xSb;cTe3_3;Se^ 
was found to present a systematic change along the insulating 



composition line (Fig. 11), which makes it possible to tune 
the properties of the surface Dirac fermions.^^^^ Such a tuning 
of the Dirac-carrier properties has been dubbed "Dirac-cone 
engineering". In particular, one can obtain both p- and ^-type 
Dirac carriers by tuning (x, y) along the insulating line, which 
is useful for designing Sip-n junction of the topological sur- 
face state. In this respect, similar band engineering has also 
been achieved in (Bi,Sb)2Te3.i^^^ 

In fact, another useful approach to reducing bulk carriers 
is to make a solid solution of Bi2Te3 and Sb2Te3. Growths 
of thin films^^^^ and nanoplates^^^^ of (Bi,Sb)2Te3 as well as 
nanoribbons^^^^ of (Bi,Sb)2Se3 have been reported, and in 
those samples the chemical potential was successfully tuned 
to be close to the Dirac point in the middle of the band gap. 
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Fig. 10. (Color online) (a) Composition-structural phase diagram of the 
Bi2-xSb;cTe3_3;Se3; system. The circles denote the compositions where the 
compensation between n- and p-type carriers are maximally achieved. The 
dashed line denotes the insulating composition previously suggested by Ter- 
amoto and Takayanagi [J. Phys. Chem. Solids 19 (1961) 124]. (b) Tempera- 
ture dependencies ofpxx for the series of BSTS samples at optimized compo- 
sitions. Note that the vertical axis is in a logarithmic scale. Taken from Ref. 
103; copyright American Physical Society (2011). 
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Fig. 11. (Color online) (a) ARPES intensities of Bi2-jcSbjcTe3_^Se^ for, 
from left to right, (x,y) = (0, 1), (0.25, 1.15), (0.5, 1.3), and (1, 2), mea- 
sured with /zy = 58 eV along the ky axis (f - ^ cut) at T = 30 K. White 
arrows and thick (red) dashed lines indicate the energy positions of the Dirac 
point (Ej)p) and the valence band top (Eyb), respectively, (b) Schematic band 
diagrams for Bi2-xSb;cTe3_3;Se3; at the insulating compositions, derived from 
the data shown in (a). Taken from Ref. 187. 



resolution ARPES experiments. ^^^^ Also, shining a circularly- 
polarized light to the surface state selectively excites electrons 
with a particular spin polarization/ ^^^ which means that elec- 
trons with a particular momentum direction is photo-excited 
in the TI surface state, and as a result, a spin-polarized pho- 
tocurrent flows in a direction dictated by the light polarization. 
Such a peculiar photocurrent has actually been experimentally 
observed. ^^^^ 

However, detection of a spin-polarized current on the sur- 
face of a TI turns out to be extremely difficult. This is because 
the spin-momentum locking makes the charge and spin scat- 
tering times to be the same, and therefore the spin diff'usion 
length is equal to the electron mean free path in the TI sur- 
face state. ^^^^ In such a situation, the spin polarization is sig- 
nificantly diminished in a diff'usive transport (to the order of 
Ak/kp, where A^ is the shift of the Fermi surface induced by 
the applied electric field), ^^^^ and the experiment must be done 
in a ballistic transport regime to detect a spin-polarized cur- 
rent. ^^^^ No such experiment has been reported for 3D TIs, but 
in the CdTe/HgTe/CdTe quantum well in the 2D TI regime, 
spin polarization of the edge current has been confirmed,^^^ 
thanks to the long electron mean free path achievable in HgTe 
quantum wells. 



8. Properties of 3D TI Materials 

8. 1 Surface state and helical spin polarization 

The most prominent property of a TI is the existence of a 
gapless surface state. The gapless nature is protected by TRS 
in Z2 topological insulators. What makes this surface state 
distinct from ordinary surface states (including accumulation 
and inversion layers) is its helical spin polarization (Fig. 1), 
which is also called spin-momentum locking; namely, the sur- 
face state is spin non-degenerate and the direction of the spin 
is perpendicular to the momentum vector and is confined in 
the surface plane. In fact, if a band has such a peculiar spin 
polarization and the system preserves TRS, there must be a 
Kramers partner for each band and Kramers theorem dictates 
that the two bands cross each other at TRIMs, which guaran- 
tees the gapless nature of the surface state [see Fig. 5(b)]. 

The helical spin polarization of the surface state means that 
a dissipationless spin current exists on the surface in equilib- 
rium, because there is not net charge flow but the spin angular 
momentum flows in the direction perpendicular to the spin di- 
rection. The helicity of the surface state (i.e. whether the up 
spin is associated with +k or -k) determines the spin current 
direction. 

The spin-momentum locking naturally gives rise to various 
interesting spin-related physics. For example, charge fluctua- 
tions are naturally accompanied by spin fluctuations, leading 
to novel spin-plasmon excitations. ^^^^ The characteristic en- 
ergy scale of such excitations are predicted to be a few meV, 
and those excitations may have been observed in ultra-high- 



8. 2 Dirac fermion physics 

When two spin-non-degenerate bands forming a Kramers 
pair cross each other at a TRIM, the energy dispersion near 
this crossing point can be approximated by a linear dispersion 
and the low energy physics is described by massless Dirac 
equation. This means that a Kramers pair of surface bands are 
actually forming a single 2D Dirac cone on which the spin 
degeneracy is lifted [Fig. 1(d)]. 

The Dirac equation for a free particle with mass m is written 
as 



Eil/{Y) = c 
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where a is the vector of Pauli matrices and p is the momen- 
tum operator. The energy eigenvalue of this equation is 



±c 



^ 



(111) 



and hence the Dirac equation always has positive and negative 
energy states. One can easily see that this energy eigenvalue 
has a gap for a finite mass m, but it becomes gapless when 
m = 0. This is the reason why a gapless system obeying the 
Dirac equation is called "massless"; it does not mean that the 
eff'ective mass of electrons becomes zero. In fact, when the 
dispersion is linear, its second derivative is zero and the eflfec- 
tive mass m* = h^(d^E/dk^)~^ diverges. 
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Graphene has emerged as a prototypical material to host 
2D Dirac fermions.^^^ Although both graphene and TIs are 
Dirac systems, there is an important difference. Namely, the 
former has both spin and valley {K and K' points in the BZ) 
degeneracies, while the latter is non-degenerate, and hence 
the Dirac fermion physics is simpler in TIs. This difference in 
degeneracy also means that the fermion degrees of freedom is 
1/4 in TIs compared to that in graphene. 

A prominent property of Dirac fermions is that they carry 
the Berry phase of n, as was initially noted by Ando et al^"^^ 
To see this, let us consider 2D massless Dirac fermions with 
the Fermi velocity vp, for which the 4 x 4 matrix equation 
reduces to 2 x 2 and the Dirac equation is written as 



Eil/{Y) = TivfCT ' kj/^(r) = -ifivpcr • Vj/^(r). 

The eigenfunctions of this equation are 

I I -ie(k)/2 \ 



(112) 



^±(r) = 



je(k)/2 



e*"" = M±(k)e*^ (113) 



V2I ±^" 
where ^(k) = arctan(^3;/^;c), and the energy eigenvalues are 

E^ = ±nvFk. (114) 

When the wave vector k is adiabatically rotated anticlockwise 
along a closed path C to encircle the origin, the Berry phase 
y acquired during this adiabatic cycle is 



y: 



i 



dk ' i{u+(k) \S/k\ u+(k)} = n. 



(115) 



Because of this n Berry phase, the time-reversed scattering 
paths, which in ordinary metals interfere constructively to 
cause the weak localization effect, ^^^^ now destructively inter- 
fere each other, leading to the weak anti-localization effect."^^^ 
Another prominent property of Dirac fermions is their pe- 
culiar Landau quantization of the energy states in magnetic 
fields. It has been shown ^^^^ that the quantization occurs in 
the manner 

E±(N) = ± ^(2envlB/c)N, (116) 

where N = 0, 1,2, ... This means that in Dirac systems when 
the Fermi level is varied, the Landau-level (LL) spacing is 
not a constant but changes as Va^, as opposed to ordinary 
metals in which the LL spacing is simply Hojc (= eHB/rricC 
where rric is the cyclotron mass) and is independent of the 
Fermi level. Also, Eq. (116) indicates that there is the zeroth 
LL with N = 0, which is pinned to the charge neutrality point 
(Dirac point). Therefore, the Landau quantization of massless 
Dirac fermions is characterized by the occurrence of the zero- 
energy state and the symmetrical appearance of Va^ states 
on both the positive and negative energy sides of the Dirac 
point [see Fig. 12(b,c)]. Such a peculiar Landau quantization 
of the surface state has indeed been observed in TIs by STS 
experiments in magnetic fields.^^'^^^'^^^^ 

When the Landau quantization of the form Eq. (116) takes 
place, the associated quantum Hall effect becomes unusual. 



and the Hall plateau between the Mh and (A^ -\- l)th LLs is 
quantized to 

^2 



xy 



e- I 1 

T^^^2 



(117) 



which is called half-integer quantization [see Fig. 12(d)]. ^^^ 
This unusual quantization can be understood to be a result 
of the existence of the N = LL at the Dirac point, which 
dictates that the first Hall plateaus on the positive and negative 
energy sides must appear antisymmetrically, because CTxy is an 
odd function of energy. The half-integer quantization can also 
be understood to be a result of the n Berry phase, ^^^'^^^^ which 
can be most easily understood by remembering the Laughlin's 
gauge argument^^^^ for the quantization of CTxy. 

8.3 Quantum oscillations 

When the Landau quantization of the energy states occurs 
in crystalline solids in magnetic fields, the density of states 
(DOS) becomes periodically modulated as a function of mag- 
netic field, which leads to various sorts of oscillation phenom- 
ena that are generally called quantum oscillations. "^^^ In par- 
ticular, the oscillations occurring in conductivity are called 
Shubnikov-de Haas (SdH) oscillations, while the oscillations 
occurring in magnetic susceptibility are called de Haas-van 
Alphen (dHvA) oscillations. 

SdH oscillations play particularly important roles in the 
studies of 3D TIs for two reasons. First, they provide a means 
to selectively and quantitatively characterize the 2D surface 
states that coexist with a 3D bulk state. Second, the phase fac- 
tor of the oscillations directly reflects the Berry phase of the 
system, which allows us to elucidate whether the electrons 
showing the SdH oscillations are Dirac fermions or not. 

In SdH oscillations, conductivity oscillates periodically as 
a function of l/B and the oscillatory part of the longitudinal 
conductivity o-xx follows 



ACTr 



COS 



2.1^.^-^ 



(118) 



where F is the oscillation frequency and /S accounts for a 
phase shift (0 < /? < 1). In fact, the same phase factor ap- 
pears in the Onsager semiclassical quantization condition"*^' 



2?re„/,, 1 

^^r"2 



(119) 



which is satisfied when the Nth LL is crossing the Fermi en- 
ergy Ef. Here, A^ is the area enclosed by electrons in the 
^- space with their cyclotron orbits on the Fermi surface. The 
parameter 13 is simply the Berry phase y divided by 2n. For 
spinless fermions, it is known"^^'^^^^ that the Berry phase is 
zero for a parabolic energy dispersion (J3 = 0) and, as already 
noted, 71 for Dirac fermions possessing a linear energy disper- 
sion (J3 = |). In real TI materials, the Dirac dispersion is not 
strictly linear but contain a parabolic component. Neverthe- 
less, it has been elucidated thatyS is robustly ^ at least at large 
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j^ 203-205) 

When Ef lies at the center of a LL (which is usually broad- 
ened due to thermal fluctuations and disorder), the DOS takes 
a maximum; on the contrary, the DOS takes a minimum when 
Ep lies in between two neighboring LLs [Fig. 12(b,c)]. In the 
latter situation, a certain number of LLs are completely filled 
and the next LL is empty. Therefore, a minimum in cTxx, which 
occurs when DOS takes a minimum, signifies a complete fill- 
ing of some N LLs, and one can assign an integer index N 
to that minimum. In ordinary metals, this LL index N cor- 
responds to the filling factor y. This can be easily understood 
by remembering the situation in the ordinary quantum Hall ef- 
fect, in which cTxy is quantized to ve^/h and cTxx becomes zero 
when the chemical potential lies between the vth and (v -i- l)th 
LLs. (In the quantum Hall eff'ect, the vanishing cTxx is a reflec- 
tion of the bulk gap opening at the Fermi level.) On the other 
hand, in the case of Dirac fermions, the filling factor y is not 
N but is N -\- ^ because of the half-integer quantization which 
essentially stems from the existence of the zeroth LL. 

The phase factor JS in the SdH oscillations can be experi- 
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Fig. 12. (Color online) (a) Partially-filled 2D Dirac cone; // is the Fermi 
level, (b) Landau quantization of the Dirac cone; LLs below // are filled with 
electrons. Note that the spacing between LLs change as Va^ and the A^ = 
LL is pinned to the Dirac point, (c) In a higher magnetic field, the spacing 
between LLs increases as ^/B, and fewer LLs are filled, (d) Schematic be- 
havior of (Txx and (Txy in the quantum Hall regime of 2D Dirac fermions; the 
situations depicted in panels (b) and (c) are marked by arrows. Notice that the 
minima in ctxx marks a complete filling of up to Mh LLs. 




mentally determined from an analysis of the so-called LL fan 
diagram, in which the sequence of the values of l/B^ corre- 
sponding to the Mh minimum in cTxx are plotted versus N. 
From Eq. (118), one can see that the Mh minimum occurs 
when the argument of the cosine equals (2N - l)7i, i.e. 






(120) 



Therefore, the plot of l/B^ vs N makes a straight line with a 
slope F corresponding to the oscillation frequency. Note that 
F/Bi = p holds for the 1st minimum. Also, it follows from 
Eq. (120) that, when a linear fit to the LL fan diagram is ex- 
trapolated to l/^A^ -^ 0, the intercept on the A/^-index axis 
gives the phase factor yS. When the yS value thus obtained is ^, 
one may conclude that the SdH oscillations come from Dirac 
fermions. 

In the past, there have confusions about the proper way 
to construct the LL fan diagram; namely, whether to assign 
an integer index to a minimum in p^x or to a minimum in 
^^^99,100,104,146,147,160,206-208) ^'^^^^ ^^ ^^ important to clar- 
ify this confusion, let us discuss this issue in some detail. In 
solids, the resistivity tensor is an inverse of the conductivity 
tensor, and in the isotropic case their relation is 

-1 



Pxx 
pyx 



Pxy 
Pxx 



O-x. 
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xy 
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xy 
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xy 
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xy 



(121) 



Therefore, when the condition cFxx ^ CTxy is satisfied (which 
is usually the case with low-carrier-density semiconductors), 
Pxx - o-xxlcy^ and the minima in pxx coincide with those 
in (J XX' This is the reason why the LL fan diagram con- 
structed from the pxx data can give the correct phase factor 
in graphene.^^^'^^^^ However, in the case of TIs, due to the 
presence of the bulk transport channel, often the condition 
o-xx ^ o-xy does not strictly hold. In the extreme case, when 
o-xx ^ o-xy (which is usually the case with metals), pxx - o-^l 
and the minima in pxx now coincide with the maxima in axx- 
Therefore, unless cFxx ^ o^xy is satisfied, performing the LL 
fan diagram analysis using the resistivity data is dangerous. 
Ideally, one should measure both Pxx(B) and Pyx(B) at the 
same time and convert them into crxx(B) and crxy(B) to per- 
form reliable LL fan diagram analysis to elucidate the correct 
Berry phase. 

In this regard, in the early stage of the TI research, the LL 
fan diagram analyses of the SdH oscillations observed in TIs 
were influenced by the case with graphene and used the min- 
ima in Pxx for indexing integer a^.99,ioo,io4,146,147,i60,206-208) 

Therefore, the conclusions regarding the Berry phase in those 
early publications should be taken with care. It was Xiong 
et al who first switched to correctly using the minima in 
(Txx for the LL fan diagram analysis, ^^^^ and some of the re- 
cent works performed reliable analyses and unambiguously 
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demonstrated the Dirac nature of the surface state through de- 
termination of the Berry phase.^^^'^^^'^^^'^^^^ In addition, to 
determine the Berry phase from the LL fan diagram in a most 
reUable way, one should fix the slope of the linear fitting by 
using the frequency F obtained from the Fourier analysis of 
the data; this way, the intercept on the A^-index axis is ob- 
tained with little ambiguity. ^^^' ^^^^ An example of the LL fan 
diagram obtained for an MBE-grown Bi2Se3 film is shown in 
Fig. 13. 

For the phase analysis of the SdH oscillations, an interest- 
ing lesson can be learned from the data for Bi2Te2Se pub- 
lished in 2010 by Ren et al^^'' In their original paper, they as- 
signed integer index N to the minima in p^^ and obtained an 
intermediate phase factor /? of 0.22 ± 0.12 [see Fig. 14(b)]. ^^^ 
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Fig. 13. (Color online) (a) SdH oscillations observed in a 10-nm- thick 
MBE-grown Bi2Se3 film; to enhance the visibility of the oscillations, 
dRyx/dB is plotted vs l/Bj_ (= l/BcosO) here. Notice that the positions of 
maxima and minima do not change with l/B±, which gives evidence for the 
2D origin of the oscillations, (b) LL fan diagram constructed from the anal- 
ysis of the (Txx data measured at 1.6 K and 6 = 0°; here, integer index N are 
assigned to the minima (upper inset). Upon making a linear fitting to the data, 
the slope is fixed at the frequency F obtained from the Fourier analysis of the 
data shown in the lower inset; the straight-line fitting extrapolates to 0.40 ± 
0.04 on the A/^-index axis, which is close to the value 1/2 expected for Dirac 
fermions. Based on the data from Ref. 149. 



However, if the data are re-analyzed by calculating ctxx from 
o'xx = P XX /(fix + Pyx) ^^^ identifying its minima to signify 
integer N, the same data set now gives fS = 0.5 [Fig. 14(c)], 
which reaffirms that Ren et al were indeed observing SdH 
oscillations from Dirac fermions. 
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Fig. 14. (Color online) (a) SdH oscillations in dpy^jdB observed in a 260- 
//m-thick Bi2Te2Se single crystal; inset shows the magnetic-field angle de- 
pendence of the oscillation frequency F . (b) LL fan diagram based on the 
dpxxIdB data at 1.6 K and ^ = 0° shown in the inset. Assuming that integer 
index N are to be assigned to the minima in pxx, the minima and maxima in 
dpxxIdB correspond ioN+\ and N+\, respectively. Linear fitting to the data 
extrapolates to 0.22 ± 0.12. (c) LL fan diagram based on the dcTxxIdB data 
shown in the inset. Linear fitting to the data (in which the slope is constrained 
from the Fourier analysis) extrapolates to 0.5. Based on the data from Ref. 99. 
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Besides the Berry phase, the SdH oscillations contain var- 
ious useful information. The Onsager's relation^^^ gives F in 
terms of the Fermi wave vector kp as 



F = (flc/27ie)7ik^, 



(122) 



and therefore one can calculate kp, and hence the carrier den- 
sity, from the frequency of the SdH oscillations. To be ex- 
plicit, for the topological 2D surface state the surface carrier 
density ris = v^^kp is directly obtained from F. If the sur- 
face state is of non-topological origin (such as accumulation 
or inversion layer), rig should be multiplied by 2 to account 
for the spin degeneracy. Also, if the SdH oscillations come 
from a 3D bulk state and the bulk Fermi surface is an ellip- 
soid, measurements of the SdH frequencies in magnetic fields 
along the three principal axes give the lengths of the three 
semiaxes, Z:^, Z:^, and k^, the bulk carrier density n^D for this 
Fermi surface is obtained as ^3d = t^y^f^f^F' where the 
spin degeneracy is taken into account. ^^^^ 

The 2D nature of the SdH oscillations can be tested by tak- 
ing the dependence of F on the angle 6 between the mag- 
netic field direction and the surface plane normal [Figs. 13(a) 
and 14(a)]; if the measured F changes as ~ 1/cos^, it is a 
strong indication that the SdH oscillations come from a 2D 
system.^5'99'100'104,146,149) Qne should note, however, that a 

reasonably wide range of 6 (e.g., up to 50°) had better be mea- 
sured, because SdH oscillations from an elongated 3D Fermi 
surface can present an approximate 1/ cos 6 scaling for a lim- 
ited range of ^.^^^^ Also, disappearance of the SdH oscillations 
for 6 = 90° gives additional support to the 2D nature. 

The cyclotron mass m^ of the carriers can be determined 
from the analysis of the temperature dependence of the am- 
plitude of the SdH oscillations. This is because the quantum 
oscillations are expressed in the Lifshitz-Kosevich theory "^^^ 
as 



Acr^j, = AqRtRdRs cos 



2.lf.i-^ 



(123) 



where Aq is a constant, Rt = 

In^ikBTinoJc)! ^mh[27i^{kBT Ifiojc)], Rd 

Qxp[-2n^(kBTD/fioJc)], and Rs = cos(^7igme/mc) are 
temperature. Dingle, and spin damping factors, respec- 
tively,"^^^ and To is called Dingle temperature (g is the 
electron g-f actor and m^ is the free electron mass). In a fixed 
magnetic field, Rd does not change and the temperature de- 
pendence show up only through Rt ; therefore, a fitting of the 
temperature dependence of the oscillation amplitude allows 
us to determine ojc, which in turn gives rric = eB/(c<jJc). 
Remember, m^ is defined as 

_ n^ i dA(E) \ 

27r\ dE 1^^^/ 



(124) 



where A(E) is the area enclosed by the cyclotron orbit in the 
/:-space. Since the cyclotron orbit is confined on the Fermi 
surface in the ^- space, the enclosed area changes with E^. 



In the case of 2D Dirac fermions with the energy dispersion 
E(k) = fiypk, one obtains A^Ep) = nkj. = 7iEp/(fiVF)^ and 
hence rric = Ef/vj^ = Mf/vf. Thus, once rric is determined 
from the temperature dependence of the SdH amplitude, one 
can calculate Vf = fikF/rric and compare it with the slope of 
the Dirac dispersion known from the ARPES measurement to 
judge if the obtained SdH result is consistent with the known 
dispersion. 

After rUc is determined, one can further determine the 
Dingle temperature T^ (= H/lnkBr) from the magnetic-field 
dependence of the amplitude of the SdH oscillations at fixed 
temperature and obtain the quantum scattering time r. This 
is usually done by making the Dingle analysis, in which 
one plots ln[AB smh(aT / B)] against l/B (A is the observed 
magnetic-field-dependent amplitude of the oscillations and 
a = 14.7 (rric /tfie) [T/K]). Such a plot makes a straight line and 
its slope gives To because of the relation A = AqRtRdRs ~ 
[(In^kBT Ifiojc)! ^mhiln^kBT Ifiojc)] exp(-27i^kBTD/ficOc)- 
The obtained r is used for calculating the mean free path 
^sdH _ ^^^^ which in turn gives an estimate of the surface 
carrier mobility yu^^^ = er/rric = e£/fikF. 

Historically, the first observation of quantum oscillations 
coming from the 2D state of a 3D TI was made by Taskin and 
Ando in 2009 for Bii_;,Sb;„ in which both dHvA and SdH 
oscillations were clearly observed.^^^ In their experiment, up 
to five frequencies were resolved in the Fourier transform of 
the oscillations, and detailed angular dependence measure- 
ments of those frequencies in all three basal planes eluci- 
dated the existence of one circular 2D Fermi surface on the 
Ci plane and three ellipsoidal electron pockets located at the 
L points in the 3D BZ. Furthermore, it was possible to deter- 
mine rric - O.OOSlrric and T^ = 6.7 K for the 2D channel from 
the analyses of the 2D component in the data. The 2D carrier 
density obtained from the oscillation frequency was only rig 
= 1.4 X 10^^ cm"^. The extremely small rric at very low ris is 
a characteristic feature of massless Dirac fermions, in which 
rric = Ef/vJ.', furthermore, the phase of the 2D component of 
the oscillations indicated the Berry phase of n,^^^'' and hence 
the observed 2D carriers were most likely of topological ori- 
gin. Nevertheless, it turned out that the 2D component of the 
oscillations is too strong if one assumes that the 2D carriers 
reside only on the outer surface of the sample, because the 
total number of bulk electrons was ~10^ times larger than the 
total number of 2D carriers, and yet, the dHvA oscillations 
from the 2D carriers were as strong as those from bulk carri- 
ers. Although the actual situation in Bii_;cSb;c is still not clear, 
what is probably happening is that Bii_;cSb;c crystals contain 
a high density of planar defects along the crystallographic Ci 
plane, and topological 2D states reside on such internal "sur- 
faces". 

The first observations of SdH oscillations coming really 
from the outer surface of 3D TIs were made for Bi2Te3 and 
Bi2Se3 by Qu et al ^^^^ and by Analytis et al}^^^ respectively, 
at nearly the same time in 2010. The 2D nature of the SdH os- 
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cillations were confirmed by taking the angular dependence of 
the oscillation frequency, as was the case with Bii_;cSb;c-^^'^^^ 
Since the surface band structures were already known from 
ARPES for those materials 7^~^^^ the agreement of vp esti- 
mated from the SdH data with that known from ARPES gave 
confidence in the origin of the observed 2D SdH oscillations. 
Those observations were made possible by minimizing the 
naturally-doped bulk carriers in Bi2Te3 and Bi2Se3, but still, 
the surface transport accounted for only less than 0.3% of the 
total conductance. 

8.4 Two-band analysis 

In real TI materials, the bulk transport channel is usually 
not negligible and one should consider parallel conductions 
through surface and bulk. Such a situation can be treated by 
a two-band model, in which the composite Hall resistivity is 
described as^^^ 



Pyx 



(RspI + RbP^JB -h RsRbiRs + Rb)B^ 



(125) 



(Ps^Pb)^^(Rs^Rb)^B^ 

where R^ and pt are the Hall coefficient and resistivity of the 
bulk channel, Rg = t/(ens), and Ps = put, with t the sample 
thickness and pn the surface sheet resistance. By fitting the 
Pyx(B) data with this two-band model, one can obtain n^D (= 
t/eRh), ph, ns, andpn, from which the mobilities for the bulk 
and surface channels are also obtained. 

One should note that this analysis involves as many as four 
parameters. Nevertheless, the parameters must be consistent 
with the pxx value in zero field, which gives one constraint. 
If SdH oscillations are observed in the surface transport, the 
SdH frequency F fixes the value of ns, which is an additional 
constraint. With those two constraints, the two-band analysis 
becomes reasonably reliable. 

For example, in the transport study of a 260-yum-thick 
BiTe2Se crystal reported by Ren et al,^^^ observation of sur- 
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Fig. 15. (Color online) Two-band analysis of the non-linear Pyx(B) data 
observed in Bi2Te2Se at 1.6 K. Open circles are the data and the solid line 
is the fitting result. In this fitting, the surface carrier density ris is not a free 
parameter but is constrained from the SdH frequency. Based on the data from 
Ref. 99. 



face SdH oscillations (Fig. 14) allowed them to fix n^, and 
the two-band analysis of pyx(B) (Fig. 15) gave consistent es- 
timates of all the transport parameters: surface mobility //^ = 
1450 cmVVs, bulk mobility ///, = 11 cm^A^s, ^, = 1.5 x 10^^ 
cm"^, and n^j) = 2.4 x 10^^ cm"-^. These values are used for 
calculating the fraction of the surface contribution in the total 
conductance, which, in this example, was ~6%. 

However, if the SdH oscillations are not observed, one 
should take the result of the two-band analysis with a 
grain of salt. This is because the model does not consider 
any magnetic-field dependence of the scattering rates, even 
though the magnetoresistancep;c;c(^) is usually large and com- 
plicated in TIs. In fact, it is almost always the case that the 
combination of parameters obtained from the two-band anal- 
yses of pyx(B) does not correctly reproduce the pxx(B) behav- 
ior. 

It is often observed in TIs that the surface mobility jj^^ es- 
timated from the two-band analysis is higher than jj,^^^ esti- 
mated from SdH oscillations. This discrepancy is most likely 
intrinsic and is a manifestation of the peculiar helical spin po- 
larization. This is because r obtained from SdH oscillations 
reflects scattering events in all directions equally, whereas in 
the transport coefficients such as pyx the backward scattering, 
which is discouraged in the TI surface states, plays the most 
important role. More specifically, r^^ to govern the resistivity 
acquires the additional factor 1/(1 - cos 0) upon spatial aver- 
aging (0 is the scattering angle), while r to govern the dephas- 
ing in the quantum oscillations is given by a simple spatial 
averaging without such a factor.^ ^^^ Hence, if the small-angle 
scattering becomes dominant (which is often the case at low 
temperature), r^^ can be much larger than r. 

8.5 Weak anti-localization effect 

As briefly mentioned in conjunction with the Dirac fermion 
physics, the n Berry phase associated with charge carriers 
leads to the weak anti-localization efl'ect.^^^ In the ordinary 
weak localization efl'ect,^^^^ electrons in metals are localized 
due to constructive interference of the electron wave func- 
tions between two time-reversed paths, which enhances the 
probability of those electrons to be found at a certain spatial 
location and reduces their ability to transport currents. This 
interference can be destroyed by breaking TRS with applied 
magnetic field, which shifts the phase factor of the wave func- 
tions of the two time-reversed paths diff'erently. Therefore, 
when the weak localization is taking place, application of a 
magnetic field leads to a negative magnetoresistance, i.e., an 
enhancement in the conductivity. The anti-localization eff'ect 
is just the opposite of this weak localization effect;^^^ namely, 
because of the phase shift of n which occurs when an elec- 
tron travels along a closed path, the two time-reversed paths 
interfere destructively, reducing the probability of electrons to 
localize. Naturally, application of a magnetic field tends to lift 
this destructive interference, leading to a positive magnetore- 
sistance. 
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It should be noted that this weak anti-locaUzation effect 
also occurs in a system with strong SOC,^^^^ which causes 
a spin rotation whenever an electron is scattered from an im- 
purity, and the resulting phase change of the electron wave 
function leads to the destructive interference similar to the 
case of the n Berry phase. Therefore, the observation of the 
weak anti-localization effect alone does not give evidence for 
the existence of Dirac fermions. 

The magnetic-field dependence of the conductivity for the 
weak anti-localization effect in 2D systems was calculated by 
Hikami, Larkin, and Nagaoka as^^^^ 
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where ^ is the digamma function and L^ is the phase coher- 
ence length. The prefactor a should be - 1 for each transport 
channel that either carries a n Berry phase^^^ or bears a strong 

In Bi2Se3 thin films, this weak anti-localization effect is 
frequently observed and a ^ -^ is usually found. This is 
because the top and bottom surfaces are connected though 
the bulk channel to form a single diffusive transport chan- 
nel. ^^^'^^^~^^^^ Remember, as long as the thickness is shorter 
than L(j) which is usually of the order of 100-1000 nm in TIs, 
electrons can diffusively travel from the bottom surface state 
to the top surface state though the bulk state using three differ- 
ent Fermi surfaces without losing the phase memory. Figure 
16 shows an example of the weak anti-localization behavior 
observed in a series of MBE-grown Bi2Se3 thin films with 
varying thickness, in which a ^ -\ was consistently observed 
unless the films are not too thin. 

Intriguingly, in films with low bulk-carrier density, by 



0.0 



< 



-0.5 



I ■__, 1 ■ ' . 1 




|\^^^^3^ 


(a) /I 


i\ 


yyl 


iV x^^s 


^ / / 


\\ ' 


'/'' 


\\ ^ 



-0.1 0.0 0.1 

B(T) 




10 100 

t(nm) 



Fig. 16. (Color online) (a) Weak anti-localization behavior observed at 1 .6 
K in (Txx of a series of MBE-grown Bi2Se3 thin films with various thickness 
shown in QL unit (1 QL is 0.95-nm thick); dashed lines are the fittings using 
the Hikami-Larkin-Nagaoka formula [Eq. (126)]. (b) Thickness dependence 
of the parameter a in Eq. (126). Inset shows schematic energy bands above 
and below the critical thickness at which hybridization of top and bottom 
surface states starts to develop. Taken from Ref. 149; copyright American 
Physical Society (2012). 



applying a gate voltage one can decouple one of the sur- 
face states from the rest through creation of an intervening 
gapped layer, and in such a situation a has been found to 
become -1, signifying two independent diffusive transport 
channels. ^^^'^^^^ 

8. 6 Topological protection of the surface state 

An important property of the topological surface state of 
a TI is its topological protection. There are three aspects in 
the notion of topological protection. One is a consequence of 
the fundamental Z2 topology which guarantees the existence 
of a gapless surface states as long as TRS is preserved. An- 
other is a consequence of the helical spin polarization, which 
makes the spin eigenvalues of the states with k and -k to 
be exactly opposite; as a result, electrons with momentum k 
cannot be backscattered into the -k state because of the spin 
mismatch [Fig. 1(d)], and hence the surface state is protected 
from backscattering. (Note, however, that the scattering prob- 
ability from k to -k -h ^ is non-zero, and hence the surface 
state is not dissipationless.) The third is the n Berry phase 
associated with massless Dirac fermions,^^^ which protects 
them from weak localization through destructive interference 
of time-reversed paths. Those three effects collectively protect 
the topological surface state. 

However, when TIs are thinned to the extent that the wave 
functions of the top and bottom surface states overlap, their 
hybridization leads to a gap opening at the Dirac point and re- 
sults in a degenerate, massive Dirac dispersion. ^^^'^^^"^^^^ For 
Bi2Se3, it was shown by ARPES experiments ^^^'^^^^ that this 
hybridization gap opens in ultrathin films with the thickness 
of 6 nm or less. 

Intriguingly, it was found by Taskin et al in transport mea- 
surements of a series of high-quality Bi2Se3 thin films that 
the opening of the hybridization gap in the Dirac dispersion 
leads to a drastic diminishment of the metallic surface trans- 
port, which helps us understand the importance of topologi- 
cal protection for the transport properties of TIs.^"^^^ Namely, 
in Fig. 16 one can see that the weak anti-localization behav- 
ior is quickly diminished in ultrathin films with thickness of 
less than 5 nm, and this is because in such ultrathin films the 
opening of the hybridization gap causes the Berry phase y to 
be reduced from n. In the simplest case, y depends on the hy- 
bridization gap A as^^^^ 



y 



"-^ 



(127) 



This change in the Berry phase diminishes the weak anti- 
localization effect and the ordinary weak localization effect 
takes over. Indeed, in Ref. 149, Pxx(T) of the 2- and 3-nm- 
thick films presented a localization behavior at low temper- 
ature, whereas all the films with thickness larger than 6 nm 
preserved a high surface mobility and presented pronounced 
SdH oscillations of Dirac nature. Therefore, the change in the 
transport properties in the hybridized regime signifies the con- 
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sequences of the loss of topological protection. 



149) 



8. 7 Magnetic Topological Insulator 

Magnetic-ion doping to 3D TIs has been studied to see the 
effect of breaking TRS. Mn-doping to Bi2Te3 was reported 
to induce bulk ferromagnetism with Curie temeperature of 
up to 12 K with 9% Mn doping.^^^^ For Bi2Se3, it was re- 
ported that Fe and Mn dopings do not induce bulk ferro- 
magnetism, but they were found to open a small gap at the 
Dirac point in the surface Dirac cone,^^^^ possibly because of 
a ferromagnetic order that develops only on the surface. Such 
a Dirac-fermion-mediated ferromagnetism was confirmed in 
Mn-doped Bi2(Se,Te)3 thin flakes, in which the chemical po- 
tential was successfully tuned into the bulk band gap so that 
the anomalous Hall effect coming only from the surface elec- 
trons can be measured to probe their ferromagnetic order. ^^"^^ 
Recently, in Cr-doped (Bi,Sb)2Te3 thin films, similar surface 
ferromagnetism was observed and, furthermore, when the 
chemical potential was tuned to be close to the Dirac point 
where a gap opens due to the ferromagnetic order, the quan- 
tum anomalous Hall effect with cr^y quantized to e^/h was 
observed at 30 mK, which signifies the appearance of a ID 
chiral edge state as a consequence of the TRS breaking in the 
2D topological surface state. ^^^^ 

9. Prospects 

9. 1 Major issues in the TI research 

There have been great improvements in the materials prop- 
erties of TIs in the past few years to make them suitable 
for fundamental research. Now there are a couple of choices 
for bulk-insulating TI materials in which the surface trans- 
port dominates over residual bulk transport at low temper- 
ature in bulk single crystals. ^^^'^^^^ High-quality thin film 
samples presenting pronounced SdH oscillations can also be 
grown, ^"^^'^^^^ which are useful for making devices to inves- 
tigate various fundamental phenomena. In future, discoveries 
of new TI materials possessing bulk band gaps of more than 
0.6 eV would be desirable for room-temperature applications, 
but existing materials such as Sn-doped Bi2Te2Se^^^^ are good 
enough for studying the topological surface transport at tem- 
peratures below -100 K. It is expected that the focus of the 
TI research will gradually move to actually realizing novel 
topological phenomena using high-quality samples that have 
become available. In this regard, there are three major themes 
for future investigations, on which I elaborate below. 

Topological magnetoelectric effects:- The topological field 
theory for Z2 TIs^^^ leads to the appearance of an unusual E -B 
term in the Lagrangian 
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where e is dielectric constant, jd is magnetic permeability, 
a = e^/hc is the fine structure constant, and ^ = or tt (mod 
In) is a topological invariant, which takes the value n in Z2 



TIs. When 6 is allowed to take any value, the second term 
(called ^-term) describes the electrodynamics of an exotic fic- 
titious particle axion. Due to the existence of the ^-term, the 
constituent equations of a TI become 

(129) 
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where D is electric induction, P is electric polarization, and 
M is magnetization. The most important consequence of this 
peculiar electromagnetism is that in TIs where 6 = n, elec- 
tric field E induces magnetization 47rM = aE and the pro- 
portionality coefficient is a universal constant a. Similarly, 
magnetic field B induces electric polarization 47rP = aB. It 
should be noted, however, that the metallic surface state elec- 
trically short-circuit the whole surface, which makes it impos- 
sible to maintain electric field or polarization necessary for 
any magnetoelectric effects; therefore, to observe the topo- 
logical magnetoelectric effect, one should gap out the sur- 
face state to avoid the short-circuiting effect.^^"^^^ This can, 
in principle, be done by depositing a ferromagnetic insula- 
tor whose magnetization is kept perpendicular to the surface. 
However, because of the practical difficulty to achieve such a 
gapping of the surface state, the topological magnetoelectric 
effect remains to be experimentally discovered. Even with- 
out opening a gap, magneto-optical properties of the surface 
state may show topological magnetoelectric effects in terms 
of Faraday rotation quantized in integer multiples of a and 
a gigantic Kerr rotation of njl^^'' Another intriguing conse- 
quence of the ^-term to be tested in future is the appearance 
of an image magnetic monopole;^^^ namely, if the surface of a 
TI is gapped out by some means and a point charge is placed 
near the surface, the response of the TI looks as if there is a 
magnetic monopole in the TI. 

Majoranafermions:- Majorana fermions are exotic charge- 
neutral particles that are their own antiparticles.^'^^^'^^^^ 
While their existence in nature as elementary particles has 
not been confirmed since its prediction in 1937, recently its 
realizations in condensed matter as quasiparticles are attract- 
ing significant attention because of their fundamental novelty 
as well as their potential for being used as a qubit of fault- 
tolerant topological quantum computing. ^^ To realize Majo- 
rana fermions in condensed matter, one must achieve particle- 
hole symmetry in a spin-non-degenerate system.^^^'^^^^ The 
former can be easily achieved in superconductors in which 
low-energy quasiparticles obey Bogoliubov-de Gennes equa- 
tion characterized by inherent particle-hole symmetry; how- 
ever, achieving a spin-non-degenerate superconductivity is 
a difficult task.^^^^ In this context, Fu and Kane realized 
that, since the topological surface state of TIs is spin non- 
degenerate, if superconductivity could be induced in the sur- 
face state by using proximity effect from an attached super- 
conductor, the resulting superconducting state harbors Ma- 
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jorana fermions.^^ There have been a number of experimen- 
tal reports to confirm the superconducting proximity eff'ect 
in the topological surface states/"^^'^^^'^^^"^^^^ but the exis- 
tence of Majorana fermions has not been elucidated. Majo- 
rana fermions may also be found in a doped TI that becomes 
a superconductor. ^-^^^ In this case, surface Dirac fermions may 
obtain superconductivity due to the proximity eff'ect from the 
bulk and become a spin-non-degenerate superconductor; an- 
other possibility is that the bulk superconductivity in a doped 
TI is itself topologically nontrivial due to strong SOC and 
its surface harbors dispersing, massless Majorana fermions. 
In this regard, the superconductor Cu;cBi2Se3,^^^^ which is a 
doped 3D TI, has attracted a lot of attention.^^^'^"^^"^"^^^ In 
fact, CuxBi2Se3 was recently found to present signatures of 
unconventional superconductivity in its point-contact spec- 
^j,^ 242) ^^^ ^^ unconventional superconductivity in this mate- 
rial is necessarily topological for symmetry reasons. ^"^^^ More 
recently, a superconducting doped TCI, Sni_;cIn;cTe, was also 
found to present similar signature of unconventional super- 
conductivity, ^"^^^ and here again, an unconventional supercon- 
ductivity is necessarily topological. ^"^^^ Those bulk supercon- 
ducting materials are promising candidates of topological su- 
perconductors to host massless Majorana fermions. 

Spintronics device applications:- The dissipationless spin 
current that exists in the topological surface state in equilib- 
rium are expected to be useful for low energy consumption 
spintronics devices. Nevertheless, it is not clear how to ex- 
ploit the dissipationless spin current for spintronics opera- 
tions. This is obviously an important theme for both theory 
and experiment. When one breaks the equilibrium and push 
current through the surface state, the resulting current is spin 
polarized; however, as already discussed in Sec. 8.1, the ex- 
pected spin polarization is extremely small in the diff'usive 
regime and the experiment should be done in the ballistic 
transport regime, ^^^^ which has been difficult in practice. If 
successful, a useful feature of the current-induced spin polar- 
ization would be that the polarization direction can be easily 
switched by tuning the chemical potential across the Dirac 
point with gating. 

9.2 New types of topologies and electron interactions 

Widening the scope of topological materials is an im- 
portant theme. So far, topological classifications of insula- 
tors and superconductors based on three discrete symmetries 
(TR, particle-hole, and chiral) have been established."^^ Re- 
cently, as mentioned in Sec. 4.5, new topological classifica- 
tions based on point-group symmetry of the crystal lattice is 
attracting significant interest, ^"^^"^"^^^ particularly after the new 
type of topological materials called topological crystalline in- 
sulators^^' ^^^^ have been experimentally discovered.^"^' ^^^' ^^^^ 
Also, whereas it was thought to be necessary for topolog- 
ical materials to have a fully-gapped energy spectrum for 
a topological invariant to be well defined,"^^ it is becoming 
possible to conceive a nontrivial topology for gapless sys- 



tems. ^^'^'^^'^'^'^^' ^"^^'^"^^^ It is an interesting time of topologi- 
cal expansion. Naturally, experimental discoveries of concrete 
materials that are nontrivial with respect to new topologies 
will continue to be crucially important. 

Another important future direction is the merger of strong 
electron correlations and topology. If a material that becomes 
both insulating and topologically nontrivial because (9/ strong 
electron correlations is discovered, such a materials would be 
called a topological Mott insulator^"^^^ and will create an en- 
tirely new and rich field of research. There are theoretical dis- 
cussions that oxides containing 5 J transition-metal elements 
might realize such a novel state of matter,^^^^ and discoveries 
along this line would be extremely interesting. 
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